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NUMERICAL  AND  ALGEBRAIC 


Gain  in  decibels  of  P2  relative  to  P\ 

G=101ogi„{P2/Pi). 


To  within  two  percent 

(27r)^^^  «  2.5;  10;  a;  20;  a;  10^. 


Euler-Mascheroni  constant^  7  =  0.57722 


Gamma  Function  r(x  +  1)  =  xr(x): 


r(l/6)  =  5.6663 

r(3/5)  = 

r(l/6)  =  4.5908 

r(2/3)  = 

r(l/4)  =  3.6256 

r{3/4)  = 

r(l/3)  =  2.6789 

r(4/5)  = 

r(2/5)  =  2.2182 

r(5/6)  = 

r(l/2)  =  1.7725  = 

r(i)  = 

1.4892 

1.3541 

1.2254 

1.1642 

1.1288 

1.0 


Binomial  Theorem  (good  for  |  i  |<  1  or  a  =  positive  integer): 


Rothe-Hagen  identity^  (good  for  all  complex  x,  j/,  2  except  when  singular): 

+  kz\  y  /y+{n-  fc)2\ 

x  +  fcz'  k  ^  y  +  {n  —  k)z  '  n  —  k  ' 

fc  =  0 

_  X  -h  y  /X  +  y  +  n2\ 
x  +  y-hn2^  n  > 

Newberger’s  summation  formula®  [good  for  \i  nonintegral,  Re  (a  +  /5)  >  — 1] 

E(— 1)  Ja-'yn{z)J0-\--yn{z)  TT  ^  ^  ^  ^  ^  ^ 

n  +  fi  - 
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VECTOR  IDENTITIES' 


Notation:  /,  5,  are  scalars;  A,  B,  etc.,  are  vectors;  T  is  a  tensor;  /  is  the  unit 
dyad. 

{1)ABxC  =  AxBC=BCxA  =  BxC.A  =  CAkB  =  CxA.B 

(2)  A  X  (B  X  C)  =  (C  X  B)  X  A  =  (A  •  C)B  -  (A  ■  B)C 

(3)  A  X  (B  X  C)  +  B  X  (C  X  A)  +  C  X  (A  X  B)  =  0 

(4)  (A  X  B)  ■  (C  X  D)  =  (A  ■  C)(B  ■  D)  -  (A  ■  D)(B  ■  C) 

(5)  (A  X  B)  X  (C  X  D)  =  (A  X  B  •  D)C  -  (A  x  B  •  C)D 

(6)  V(/s)  =  V(9/)  =  /V9  +  9V/ 

(7)  V  ■  (/A)  =  /V  ■  A  +  A  •  V/ 

(8)  V  X  (/A)  =  /V  X  A  +  V/  X  A 

(9) V.(AXB)  =  B  VXA-A.VXB 

(10)  V  X  (A  X  B)  =  A(V  ■  B)  -  B(V  •  A)  +  (B  •  V)A  -  (A  •  V)B 

(11)  A  X  (V  X  B)  =  (VB)  •  A  -  (A  •  V)B 

(12)  V(A  ■  B)  =  A  X  (V  X  B)  +  B  X  (V  X  A)  +  (A  •  V)B  +  (B  ■  V)A 

(13)  V^f  =  V-  V/ 

(14)  V^A  =  V(V  .  A)  -  V  X  V  X  A 

(15)  V  X  V/  =  0 
(10)  V  •  V  X  A  =  0 

If  ei,  62,  63  are  orthonormal  unit  vectors,  a  second-order  tensor  T  can  be 
written  in  the  dyadic  form 

(17)  T  = 

In  cartesian  coordinates  the  divergence  of  a  tensor  is  a  vector  with  components 

(18)  (V.T)i  =  X^.OTji/Sxj) 

(This  definition  is  required  for  consistency  with  Eq.  (29)].  In  general 

(19)  V  ■  (AB)  =  (V  •  A)B  -I-  (A  •  V)B 

(20)  V  •  (/T)  =  V/-T+/V-T 
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Let  r  =  ix  +  jy  +  kz  be  the  radius  vector  of  magnitude  r,  from  the  origin  to 
the  point  x,  y,  z.  Then 

(21)  V  ■  r  =  3 

(22)  V  X  r  =  0 

(23)  Vt-  =  r/r 

(24)  V(l/r-)  =  -r/r^ 

(25)  V  ■  (r/r^)  =  47r5(r) 

(26)  Vr  =  / 

If  V  is  a  volume  enclosed  by  a  surface  S  and  dS  —  ndS,  where  n  is  the  unit 
normal  outward  from  V, 


(27) 

/  dVVf  =  /  dSf 

Jv  Js 

(28) 

/  dVV  ■  A  =  /  dS  4 
J  V  J  s 

(29) 

/  dVV  T-  dS  T 
J  V  Js 

(30) 

/  dVV  X  A  =  /  dS) 

J  V  J  s 

(31) 

/  iiV(/V=s-sVV)  = 

/  V 

(32) 

f  dV{A  ■  V  X  V  X  B  - 

-I- 

Js 


(fVg  -  gVf) 

V  X  A) 


dS-(BxVxA  —  AxVxB) 


If  S  is  an  open  surface  bounded  by  the  contour  C,  of  which  the  line  element  is 
d\, 


(33) 


dS  X  V/  = 


dlf 
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(34) 

(35) 

(36) 


J  dS  V  X  A  =  ^  dl  -  A 
j  (dS  X 

I 


V)  X  A  =  (p  d\  X  A 


dS  ■  (V/  xVg)=  Cj)  fdg  = 


gdf 


DIFFERENTIAL  OPERATORS  IN 
CURVILINEAR  COORDINATES^ 


Cylindrical  Coordinates 


Divergence 


V  •  A  = 


-  —  (rAr)  +  + 

r  dr  r  dip 


dA, 

dz 


Gradient 


(V/)r 


il. 

dr  ’ 


= 


r  dip  ’ 


(V/)* 


dz 


Curl 

(V  X  A), 


1  aAr  _  djU 

r  dp  dz 


o/ir 

(V  X  A)^  =  — 


dAz 

dr 


(V  X  A)z  = 


-§-{rA^) 


1  dAr 
r  dp 


Laplacian 


vV  =  ■ 


1 

Sr 


1  d^f  d^f 

r2 


6 


Laplacian  of  a  vector 


(V^A),  =  - 

(V^A)^  =  V=A*  + 
(V^A)^  = 


2  3A0 
d(f> 


2  dAr 
dip 


r2 


Components  of  (A  •  V)B 


(A.VB).  =  A,^  +  :^:^  +  A. 


k  dBr 

dip 


dBr 
'  dz 


T 


(A  ■  VB)*  = 

OT  T  dip  dz  r 


dr  r  dip  dz 


Divergence  of  a  tensor 

( V  .  T).  =  i  |- (rT..)  +  i  +  ^  -  5^^: 

r  dr  r  dip  dz  r 


(T7  T'\  —  ^  ^  (r-T  ^  4-  ^  9T,p^  dTztj)  ^01 

(V  •  /  )0  —  ~-^{rTr^)  + - -  +  — - - i - 

r  dr  r  dip  dz  r 


(V  .  T)z  =  --firTr-.)  +  -  ^  +  ^ 
r  dr  r  dtp  dz 
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Spherical  Coordinates 


Divergence 


V  •  A  = 


dr 


Ar)  + 


1  a 

rsinfl  86 


(sin^Afl)  + 


1  dA^ 
rsin0  d(}> 


Gradient 


(V/),.  = 


dr  ' 


i^f)o  = 


i£/. 

r  86' 


(V/)^  = 


1  df 

r  sin  6  d(f> 


Curl 

(V  X  A), 


1  d 
r  sin  6  dO 


(sin  6A^) 


1  dAo 
rsin^  d4> 


(V  X  A),,  = 


rsin^  d<f>  r  dr^ 


(V><A),  =  i|:(rA»)-i^ 


Laplacian 

r^sin^e  8.(i2 


Laplacian  of  a  vector 

«  5  2 A,.  2  dAe  2cot^A0  2  dA^ 

(V’A)r  =  V=A,. - 5^  -  - 5 - 

'  dO  r^sinp  dq> 


vV 


=  li.  (r^^\  . 

5r  V  dr  } 


1 


a 


sin  6  d$ 


ty  o  2  dAj*  Aq 

(V  A)o  =  V  i40  +  - 2  •  2  fl 

80  sm**  6 


2cos0  dA^ 
r*  sin^  6  d<p 


dAr  2  cos  9  dAo 


^  sin  0  dtp  r^  sin^  6 
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Components  of  (A  •  V)B 

(A  ■  VBl  =A  I  I  AeBp-i-A^B^ 

''  ^  dr  r  dd  rsin6>  d4>  r 


(A  VB')  =  A  I  _  cot  OA^B^ 

^  ^  dr  r  dO  rsin0  d(f)  r  r 


(A  VB)  =A  I  dB^,  A^  dB^  A^Br  cotgA^Bg 

^  dr  r  de  rsin9  d<f>  r  r 


Divergence  of  a  tensor 

(V  T)r=^^{r^T,.)  + 

(V  .  7)«  =  ll(r^T^,)  + 

(^■Th  =  ^^yT.,)  + 


1  d 
rsin  $  d6 

+ 

1  d 

rsinO  dO 

+ 

1  d 
rsin  9  dO 


(sin  BTffr) 

1  dT^r 
rsin  6  d4> 

(sin  BT$$) 

1  ar^e 

rsin  9  d4> 

(sin  BTo^) 


+ 


T$g  + 
r 


T$r  _  cot  9T^^ 
r  r 


+  ^  _|_  ^0j-  9T^e 

rsin0  0^  r  r 
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DIMENSIONS  AND  UNITS 

To  get  the  value  of  a  quantity  in  Gaussian  units,  multiply  the  value  ex¬ 
pressed  in  SI  units  by  the  conversion  factor.  Multiples  of  3  in  the  conversion 
factors  result  from  approximating  the  speed  of  light  c  =  2.9979  x  cm/sec 
3  X  10^*^  cm/sec. 


Dimensions 

Physical 

Sym- 

SI 

Conversion 

Quantity 

bol 

SI 

Gaussian 

Units 

Factor 

Units 

Capacitance 

C 

i  q 

ml^ 

t 

farad 

9  X  lO” 

cm 

Charge 

9 

9 

„,l/2,3/2 

coulomb 

3  X  10” 

statcoulomb 

Charge 

density 

P 

9 

P 

P/H 

coulomb 

/m^ 

3  X  10” 

statcoulomb 

/cm^ 

Conductance 

IS— 

t 

siemens 

9  X  10“ 

cm/sec 

Conductivity 

a 

iil 

1 

t 

siemens 

/m 

9  X  10” 

sec"* 

Current 

/,z 

9 

7 

ampere 

3  X  10” 

statampere 

Current 

density 

j.j 

9 

I'^t 

-^1/2 
m  ' 

ampere 

/m= 

3  X  10” 

statampere 

/cm^ 

Density 

P 

m 

m 

¥ 

kg/m^ 

10”” 

g/cm” 

Displacement 

D 

9 

m'/" 

P/H 

coulomb 

/u? 

127r  X  10” 

statcoulomb 

/cm^ 

Electric  field 

E 

ml 

^1/2 
m  ' 

v-m 

volt/m 

-  X  10"* 

3 

statvolt/cm 

Elcctro- 

motance 

Emf 

ml"^ 

t 

volt 

-  X  10“^ 

3 

statvolt 

Energy 

V,  W 

ml^ 

t2 

mP 

P 

joule 

erg 

Energy 

density 

w,  t 

m 

m 

w 

joule/m^ 

10 

erg/cm^ 
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Dimensions 


Physical 

Quantity 

Force 

Frequency 

Impedance 

Inductance 

Length 

Magnetic 

intensity 

Magnetic  flux 

Magnetic 

induction 

Magnetic 

moment 

Magnetization 

Magneto- 

motance 

Mass 

Momentum 

Momentum 

density 

Permeability 


Sym- 

SI 

Conversion 

Gaussian 

bol 

SI 

Gaussian 

Units 

Factor 

Units 

dyne 

F 

(2 

newton 

10® 

1 

1 

hertz 

1 

hertz 

t 

7 

t 

1  11 

sec/cm 

Z 

ohm 

-  X  10  ” 

tq^ 

1 

9 

L 

ml^ 

'V 

1 

henry 

ixio-” 

9 

sec^/cm 

1 

1 

1 

meter  (m) 

10^ 

centimeter 

(cm) 

H 

n 

It 

«,l/2 
m  ' 

TJh 

ampere- 

turn/m 

47r  X  10"® 

oersted 

maxwell 

weber 

10® 

tq 

^1/2 

B 

tq 

I'/^t 

tesla 

10'* 

gauss 

I'^q 

oersted- 

m, 

amperc-m^ 

10® 

cm® 

M 

Q 

It 

^1/2 
m  ' 

TTh 

ampere- 

turn/m 

10“® 

oersted 

M, 

g 

ampere- 

477 

gilbert 

t 

10 

Mmf 

turn 

771,  M 

m 

m 

kilogram 

10® 

gram  (g) 

(kg) 

kg-m/s 

g-cm/sec 

p,P 

t 

t 

10® 

m 

Jh 

m 

1h 

kg/m^-s 

10"^ 

g/cm^-sec 

henry/m 

1  7 

9= 

1 

-  X  10^ 

477 
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Dimensions 

Physical 

Sym- 

SI 

Conversion 

Gaussian 

Quantity 

bol 

SI 

Gaussian 

Units 

Factor 

Units 

Permittivity 

e 

ml^ 

1 

farad/m 

3671-  X  10® 

— 

g 

coulomb/m^ 

3  X  10® 

p 

t^/H 

/cm^ 

Potential 

V,<f> 

ml^ 

volt 

i  X  10“® 

statvolt 

t'^q 

t 

3 

Power 

P 

ml^ 

(3 

ml^ 

(3 

watt 

10^ 

erg/sec 

Power 

m 

m 

1(3 

watt/m® 

10 

erg/cm^-sec 

density 

Pressure 

p,  P 

m 

1(2 

m 

1(2 

pascal 

10 

dyne/cm^ 

Reluctance 

n 

1 

ampere-turn 

477  X  10"® 

cm 

ml'^ 

7 

/weber 

Resistance 

R 

ml^ 

( 

7 

ohm 

-  X  10-" 

9 

sec/cm 

Resistivity 

v,  P 

ml^ 

(g2 

( 

ohm-m 

i  X  10-® 

9 

sec 

Thermal  con- 

«,  k 

m/ 

ml 

watt/m- 

10® 

erg/cm-sec- 

ductivity 

(3 

(3 

deg  (K) 

deg  (K) 

t 

( 

( 

second  (s) 

1 

second  (sec) 

m/ 

weber/m 

Vector 

A 

10® 

gauss-cm 

potential 

(g 

Velocity 

V 

( 

t 

m/s 

10^ 

cm/sec 

Viscosity 

m 

Tf 

m 

77 

kg/m-s 

10 

poise 

1 

1 

Vorticity 

c 

t 

i 

s  ' 

1 

sec 

joule 

Work 

W 

(2 

<2 

lO’’ 

erg 
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INTERNATIONAL  SYSTEM  (SI)  NOMENCLATURE® 


Physical 

Name 

Symbol 

Physical 

Name 

Symbol 

Quantity 

of  Unit 

for  Unit 

Quantity 

of  Unit 

for  Unit 

*length 

meter 

m 

electric 

potential 

volt 

V 

*mass 

kilogram 

kg 

electric 

ohm 

n 

*time 

second 

s 

resistance 

*current 

ampere 

A 

electric 

conductance 

siemens 

S 

*temperature 

kelvin 

K 

electric 

farad 

F 

^amount  of 
substance 

mole 

mol 

capacitance 

magnetic  flux 

weber 

Wb 

*luminous 

candela 

cd 

intensity 

magnetic 

inductance 

henry 

H 

fplane  angle 

radian 

rad 

magnetic 

tesla 

T 

fsolid  angle 

steradian 

sr 

intensity 

frequency 

hertz 

Hz 

luminous  flux 

lumen 

Im 

energy 

joule 

J 

illuminance 

lux 

lx 

force 

newton 

N 

activity  (of  a 
radioactive 

becquerel 

Bq 

pressure 

pascal 

Pa 

source) 

power 

watt 

W 

absorbed  dose 
(of  ionizing 

gray 

Gy 

electric  charge 

coulomb 

c 

radiation) 

*SI  base  unit  fSuppIementary  unit 


METRIC  PREFIXES 


Multiple 

Prefix 

Symbol 

Multiple 

Prefix 

Symbol 

10"^ 

deci 

d 

10 

deca 

da 

10“^ 

centi 

c 

10^ 

hecto 

h 

10“® 

milH 

m 

10® 

kilo 

k 

10"® 

micro 

10® 

mega 

M 

10"® 

nano 

n 

10® 

giga 

G 

10-‘® 

pico 

P 

10*® 

tera 

T 

10"*® 

femto 

f 

10*® 

peta 

P 

10“*® 

atto 

a 

10*® 

exa 

E 
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PHYSICAL  CONSTANTS  (SI)’’ 


Physical  Quantity 

Symbol 

Value 

Units 

Boltzmann  constant 

k 

1.3807  X  10“’’’® 

JK“' 

Elementary  charge 

e 

1.6022  X  10"'° 

c 

Electron  mass 

me 

9.1094  X  10"’*‘ 

kg 

Proton  mass 

m,> 

1.6726  X  10“”’ 

kg 

Gravitational  constant 

G 

6.6726  X  10“” 

m'^s“‘'kg"‘ 

Planck  constant 

h 

6.6261  X  10“”'' 

J  s 

< 

1.0546  X  10“”'* 

Js 

Speed  of  light  in  vacuum 

c 

2.9979  X  10® 

m  s"* 

Permittivity  of 

€0 

8.8542  X  10“'” 

Fm"' 

free  space 

Permeability  of 

MO 

477  X  10“’ 

Hm“' 

free  space 

Proton/electron  mass 

TTXjy  j nXe 

1.8362  X  10® 

ratio 

Electron  charge/mass 

c/mc 

1.7588  X  10" 

Ckg“' 

ratio 

Rydberg  constant 

r. 

1.0974  X  10’ 

Bohr  radius 

ao  =  €oh‘^/nm€^ 

5.2918  X  10“" 

m 

Atomic  cross  section 

■nao^ 

8.7974  X  10“”' 

....,2 

m 

Classical  electron  radius 

re  =  /Aneomc^ 

2.8179  X  10“'® 

m 

Thomson  cross  section 

(8^/3)re’' 

6.6525  X  10“”® 

m” 

Compton  wavelength  of 

hfmeC 

2.4263  X  10“'” 

m 

electron 

hImeC 

3.8616  X  10“'” 

m 

Fine-structure  constant 

a  =  f2€ohc 

Q-^ 

7.2974  X  10"® 
137.04 

First  radiation  constant 

Cj  =  2TThc^ 

3.7418  X  10“” 

Wm” 

Second  radiation 

C2  =  hc/k 

1.4388  X  10“” 

m  K 

constant 

Stefan-Boltzmann 

a 

5.6705  X  10“® 

Wm“”K“‘‘ 

constant 
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Physical  Quantity 

Symbol 

Value 

Units 

Wavelength  associated 

Ao  =  /ic/e 

1.2398  X  10“® 

m 

with  1  eV 

Frequency  associated 

t/Q  =  e/Zi 

2.4180  X  10“ 

Hz 

with  1  eV 

Wave  number  associated 

fco  =  e//ic 

8.0655  X  10® 

m”^ 

with  1  eV 

Energy  associated  with 

huo 

1.6022  X  10““ 

J 

1  eV 

Energy  associated  with 

he 

1.9864  X  10““® 

J 

1  m~* 

Energy  associated  with 

me^ /Seo^h^ 

13.606 

eV 

1  Rydberg 

Energy  associated  with 

kje 

8.6174  X  10“® 

eV 

1  Kelvin 

Temperature  associated 

efk 

1.1604  X  10'* 

K 

with  1  eV 

Avogadro  number 

Na 

6.0221  X  10“® 

mol“^ 

Faraday  constant 

F  =  NAe 

9.6485  X  10“ 

CrnoP^ 

Gas  constant 

NAk 

8.3145 

JK-^mor‘ 

Loschmidt’s  number 

no 

2.6868  X  10“® 

(no.  density  at  STP) 

Atomic  mass  unit 

TUu 

1.6605  X  10“““ 

kg 

Standard  temperature 

To 

273.15 

K 

Atmospheric  pressure 

po  —  nokTo 

1.0133  X  10® 

Pa 

Pressure  of  1  mm  Hg 

1.3332  X  10“ 

Pa 

(1  torr) 

Molar  volume  at  STP 

Vo  —  RTq/po 

2.2414  X  10““ 

m^ 

Molar  weight  of  air 

Mair 

2.8971  X  10““ 

kg 

calorie  (cal) 

4.1868 

J 

Gravitational 

9 

9.8067 

m 

acceleration 
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Boltzmann  constant 
Elementary  charge 


Electron  mass 
Proton  mass 
Gravitational  constant 
Planck  constant 

Speed  of  light  in  vacuum 
Proton/electron  mass 


h 

h  =  h/^-n 


6.672C  X  10”® 
6.6261  X 
1.0546  X  10“^^ 
2.9979  X  10“" 
1.8362  X  10^ 


erg/deg  (K) 
statcoulomb 
(statcoul) 


dyne-cm^/g^ 


Electron  charge/mass 
ratio 

Rydberg  constant 


Atomic  cross  section 
Classical  electron  radius 
Thomson  cross  section 
Compton  wavelength  of 
electron 

Fine-structure  constant 


First  radiation  constant 
Second  radiation 
constant 

Stefan-Boltzmann 


Wavelength  associated 
with  1  eV 


statcoul/g 


-  c/l3 
ao  =  h^fme^ 
TTOO^ 

r*  =  /mc^ 

(87r/3)r,2 

h/nicC 

h/nXeC 

Q  =  f  he 
a-’ 

Cl  =  27rhc^ 

C2  =  hc/k 


1.0974  X  10® 


6.6525  X  10“''®  cm-* 

2.4263  X  10“^°  cm 
3.8616  X  10'“  cm 
7.2974  X  10”® 

137.04 

3.7418  X  10"®  erg-cm^/sec 
1.4388  cm-deg(K) 

5.6705  X  10”®  erg/cm^- 
sec-deg^ 

1.2398  X  10”*^  cm 
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Physical  Quantity 

Symbol 

Value 

Units 

Frequency  associated 

t'o 

2.4180  X  10“* 

Hz 

with  1  eV 

Wave  number  associated 

ko 

8.0655  X  10® 

cm“^ 

with  1  eV 

Energy  associated  with 

1.6022  X  10-*® 

erg 

1  eV 

Energy  associated  with 

1.9864  X  lO-*** 

erg 

1  cm”^ 

Energy  associated  with 

13.606 

eV 

1  Rydberg 

Energy  associated  with 

8.6174  X  10-® 

eV 

1  deg  Kelvin 

Temperature  associated 

1.1604  X  10'* 

deg  (K) 

with  1  eV 

Avogadro  number 

Na 

6.0221  X  10®® 

mol“^ 

Faraday  constant 

F  =  NAe 

2.8925  X  10*'* 

statcoul/mol 

Gas  constant 

R  =  NAk 

8.3145  X  10® 

erg/deg-mol 

Loschmidt’s  number 

no 

2.6868  X  10*® 

cm“® 

(no.  density  at  STP) 

Atomic  mass  unit 

THy 

1.6605  X  10"®'* 

g 

Standard  temperature 

To 

273.15 

deg{K) 

Atmospheric  pressure 

Po  =  nokTo 

1.0133  X  10® 

dyne/cm^ 

Pressure  of  1  mm  Hg 

1.3332  X  10® 

dyne/cm^ 

(1  torr) 

Molar  volume  at  STP 

Vo  =  RTq/pq 

2.2414  X  10^ 

cm^ 

Molar  weight  of  air 

Mair 

28.971 

g 

calorie  (cal) 

4.1868  X  10® 

erg 

Gravitational 

9 

980.67 

cm/sec^ 

acceleration 
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FORMULA  CONVERSION® 

Here  a.  =  10^cmm“^,  /9  —  lO^ergJ  co  —  8.8542  x  10  Fm  , 
tio  =  4TrxlO"‘'^  Hm"^  c  =  (coMo)"^^^  =  2.9979  x  10®  m  and  h  =  1.0546  x 
10“^‘‘  J  s.  To  derive  a  dimensionally  correct  SI  formula  from  one  expressed  in 
Gaussian  units,  substitute  for  each  quantity  according  to  Q  =  kQ,  where  k  is 
the  coefficient  in  the  second  column  of  the  table  corresponding  to  Q  (overbars 
denote  variables  expressed  in  Gaussian  units).  Thus,  the  formula  ao  =  h  /me 
for  the  Bohr  radius  becomes  oao  =  (h^)^/[(m.^/Q^)(e^Q:/3/47reo)]>  or  “o  = 
€oh'^/7rm.c^.  To  go  from  SI  to  natural  units  in  which  h  =  c  =  1  (distinguished 
by  a  circumflex),  use  Q  =  where  k  is  the  coefficient  corresponding  to 

Q  in  the  third  column.  Thus  oq  =  47rcoh^/[{ihh./c)(e^£ohc)]  =  Anfme  .  (In 
transforming  from  SI  units,  do  not  substitute  for  cq,  |J■o^  or  c.) 


Physical  Quantity 

Gaussian  Units  to  SI 

Natural  Units  to  SI 

Capacitance 

a/47reo 

Co“’ 

Charge 

(a/3/4^co)’^^ 

(eohe)-'/^ 

Charge  density 

(/3/47ra5eo)'^^ 

(eohe)-'/^ 

Current 

{a^/4n€o)^^^ 

Current  density 

(w/fic)*'" 

Electric  field 

(47r/3€o/a^)*^^ 

(eo/hc)^'^ 

Electric  potential 

(4ir;8co/Q)'^^ 

Electric  conductivity 

(4ir£o)"‘ 

Energy 

0 

(he)-' 

Energy  density 

/3/a^ 

(M'‘ 

Force 

01a 

Frequency 

1 

C  ^ 

Inductance 

Aitco/oc 

Length 

a 

1 

Magnetic  induction 

(4:i-/3/aVo)‘^“ 

Magnetic  intensity 

Mass 

0/a'^ 

c/h 

Momentum 

01a 

h-^ 

Power 

0 

Pressure 

0/a^ 

(Re)-* 

Resistance 

47rco/ct 

Time 

1 

c 

Velocity 

a 

c"  ^ 
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MAXWELL’S  EQUATIONS 


Name  or  Description 

SI 

Gaussian 

Faraday’s  law 

V  X  E  = - 

dt 

r,  ^  I 

V  X  E  = - — 

c  dt 

Ampere’s  law 

dD 

VXH=-+J 

„  „  1  8D  47r  , 

VxH=-  —  +  — J 
c  ot  c 

Poisson  equation 

V  D  =  p 

V  •  D  =  47rp 

[Absence  of  magnetic 
monopoles] 

VB  =0 

V  ■  B  =  0 

Lorentz  force  on 

g(E  +  v  X  B) 

9  (e+  iv  X  B  ) 

charge  q 

Constitutive 

D  =  cE 

D  =  eE 

relations 

B  =  pH 

B  =  pH 

In  a  plasma,  /i  as  =  47r  x  10  ^  H m  ^  (Gaussian  units:  ^  «  1).  The 
permittivity  satisfies  €  as  €o  =  8.8542  x  10“*^  Fm"*^  (Gaussian:  €  ft;  1) 
provided  that  all  charge  is  regarded  as  free.  Using  the  drift  approximation 
vx  =  E  X  B/B^  to  calculate  polarization  charge  density  gives  rise  to  a  dielec¬ 
tric  constant /C  =  e/eo  =  1  +  367r  x  10®p/B^  (SI)  =  H-47rpc^/B®  (Gaussian), 
where  p  is  the  mass  density. 

The  electromagnetic  energy  in  volume  V  is  given  by 


Stt  I 
Jv 


dV(H  •  B  +  E  D) 
dl^(H  •  B  +  E  D) 


Poynting’s  theorem  is 


dW 

~W 


N-dS 


L 


dVJ  •  E, 


(SI) 


(Gaussian). 


where  S  is  the  closed  surface  bounding  V  and  the  Poynting  vector  (energy  flux 
across  S)  is  given  by  N  =  E  x  H  (SI)  or  N  =  cE  X  H/47r  (Gaussian). 
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ELECTRICITY  AND  MAGNETISM 

In  the  following,  e  =  dielectric  permittivity,  ^  =  permeability  of  conduc¬ 
tor,  /i'  =  permeability  of  surrounding  medium,  a  =  conductivity,  /  =  = 

radiation  frequency,  Km  =  *^nd  =  e/eo-  Where  subscripts  are  used, 

‘1’  denotes  a  conducting  medium  and  ‘2’  a  propagating  (lossless  dielectric) 
medium,  All  units  are  SI  unless  otherwise  specified. 


Permittivity  of  free  space 
Permeability  of  free  space 

Resistance  of  free  space 

Capacity  of  parallel  plates  of  area 
A,  separated  by  distance  d 
Capacity  of  concentric  cylinders 
of  length  I,  radii  a,  6 
Capacity  of  concentric  spheres  of 
radii  a,  b 

Sclf-inductance  of  wire  of  length 
I,  carrying  uniform  current 

Mutual  inductance  of  parallel  wires 
of  length  I,  radius  a,  separated 
by  distance  d 

Inductance  of  circular  loop  of  radius 
6,  made  of  wire  of  radius  o, 
carrying  uniform  current 
Relaxation  time  in  a  lossy  medium 
Skin  depth  in  a  lossy  medium 
Wave  impedance  in  a  lossy  medium 

Transmission  coefficient  at 
conducting  surface® 

(good  only  for  T  1) 

Field  at  distance  r  from  straight  wire 
carrying  current  I  (amperes) 

Field  at  distance  z  along  axis  from 
circular  loop  of  radius  a 
carrying  current  I 


Co  =  8.8542  X  10"^^  F  m“^ 

/lo  =  ^TT  X 

=  1.2566  X  10“®Hm"^ 

Ro=  =  376.73  n 

C  =  eA/d 

C  =  2nell  ln(6/a) 

C  =  47rea6/(b  —  a) 

L  = 

L  —  {tx'l/i-n)  [1-1-4  ln(d/o)] 

L  =  -  2]  +  m/4} 

T  =  C/CT 

z  =  1m/(£  + 

T  =  4.22  X 

Bg  =  ^//27rr  tesla 

=  0.2//rgauss  (r  in  cm) 

B^  =  ,ia*//[2(a“  +2^)^/^) 
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ELECTROMAGNETIC  FREQUENCY/ 
WAVELENGTH  BANDS^° 


Designation 

Frequency  Range 

Wavelength  Range 

Lower 

Upper 

Lower 

Upper 

ULF* 

30  Hz 

10  Mm 

VF* 

30  Hz 

300  Hz 

1  Mm 

10  Mm 

ELF 

300  Hz 

3  kHz 

100  km 

1  Mm 

VLF 

3  kHz 

30  kHz 

10  km 

100  km 

LF 

30  kHz 

300  kHz 

1  km 

10  km 

MF 

300  kHz 

3  MHz 

100  m 

1  km 

HF 

3  MHz 

30  MHz 

10  m 

100  m 

VHF 

30  MHz 

300  MHz 

1  m 

10  m 

UHF 

300  MHz 

3  GHz 

10  cm 

1  m 

SHFt 

3  GHz 

30  GHz 

1  cm 

10  cm 

S 

2.6 

3.95 

7.6 

11.5 

G 

3.95 

5.85 

5.1 

7.6 

J 

5.3 

8.2 

3.7 

5.7 

H 

7.05 

10.0 

3.0 

4.25 

X 

8.2 

12.4 

2.4 

3.7 

M 

10.0 

15.0 

2.0 

3.0 

P 

12.4 

18.0 

1.67 

2.4 

K 

18.0 

26.5 

1.1 

1.67 

R 

26.6 

40.0 

0.75 

1.1 

EHF 

30  GHz 

300  GHz 

1  mm 

1  cm 

Submillimeter 

300  GHz 

3  THz 

100  /im 

1  mm 

Infrared 

3THz 

430  THz 

700  nm 

100  ^m 

Visible 

430  THz 

750  THz 

400  nm 

700  nm 

Ultraviolet 

750  THz 

30  PHz 

10  nm 

400  nm 

X  Ray 
Gamma  Ray 

30PHz 

3EHz 

3EHz 

100  pm 

10  nm 

100  pm 

In  spectroscopy  the  angstrom  is  sometimes  used  (lA  =  10  ®  cm  =  0.1  nm). 
*The  boundary  between  ULF  and  VF  (voice  frequencies)  is  variously  defined. 
fThe  SHF  (microwave)  band  is  further  subdivided  approximately  as  shown.^^ 
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AC  CIRCUITS 

For  a  resistance  R,  inductance  L,  and  capacitance  C  in  series  with 
a  voltage  source  V  =  exp(iwt)  (here  i  =  \/— i)i  the  current  is  given 
by  /  =  dqfdt,  where  q  satisfies 


^^dt^  C 


Solutions  are  q{t)  =  9a  4-  9t>  -^(0  =  where  the  steady  state  is 

la  =  iu)qa  =  V/Z  in  terms  of  the  impedance  Z  =  R-\-  i(wL  —  l/wC)  and 
=  dqtfdt.  For  initial  conditions  9(0)  =  go  =  9o  +  9s*  ^(0)  =  ^o*  the 
transients  can  be  of  three  types,  depending  on  A  =  H  —  4L/C: 

(a)  Overdamped,  A  >  0 

,,  =  +  ^  exp(-7-t)  -  exp(-7+t), 

7+  -  7-  -  7- 

^  7h-(/0  + 7-^-0)  _  ■V-(^o+7-*g-ol 

7+  -  7-  -  7- 


where  7±  =  (ft  ±  A^^^)/2L; 

(b)  Critically  damped,  A  =  0 


9t  =  [90  +  (h  +  7fl9o)tl  exp(-7fit), 

It  =  [/o  -  (h  +  7K9o)7Rtl exp(-7Rt), 


where  -^n.  =  RI‘2L\ 

(c)  Underdamped,  A  <  0 

7R90  +  ^0  ^  gocosojit  exp(-7flt), 

U),  J 

sin(c.,t)l  exp(-7Rt). 

U»1  J 

Here  wi  =  u;o(l  -  where  wq  =  (LC)"^^^  is  the  resonant 

frequency.  At  w  =  wq,  Z  =  R.  The  quality  of  the  circuit  is  Q  =  loqL/R. 
Instability  results  when  L,  R,  C  are  not  all  of  the  same  sign. 


li  =  /o  cos  Wit  - 
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DIMENSIONLESS  NUMBERS  OF  FLUID  MECHANICS*^ 


Name(s) 

Symbol 

Definition 

Significance 

Alfven, 

Karman 

Al,  Ka 

Va/V 

*(Magnetic  force/ 
inertial  force)^'^^ 

Bond 

Bd 

(p'  -  p)L^g/S 

Gravitational  force/ 
surface  tension 

Boussinesq 

B 

V/{2gRfl-^ 

(Inertial  force/ 

gravitational  force)^/^ 

Brinkman 

Br 

pV^/kAT 

Viscous  heat/conducted  heat 

Capillary 

Cp 

mV/e 

Viscous  force/surface  tension 

Carnot 

Ca 

(J2  -  Ti)/T2 

Theoretical  Carnot  cycle 
efficiency 

Cauchy, 

Hooke 

Cy.  Hk 

pV^ir  = 

Inertial  force/ 

compressibility  force 

Chandra¬ 

sekhar 

Ch 

B'^L^/pvg 

Magnetic  force/dissipative 
forces 

Clausius 

Cl 

LV^p/kAT 

Kinetic  energy  flow  rate/heat 
conduction  rate 

Cowling 

C 

(VA/Vf  =  Al^ 

Magnetic  force/inertial  force 

Crispation 

Cr 

Effect  of  diffusion/effect  of 
surface  tension 

Dean 

D 

Transverse  flow  due  to 

curvature/longitudinal  flow 

[Drag 

coefficient] 

Cd 

(p'  -  p)Lg/ 
p'V^ 

Drag  force/inertial  force 

Eckert 

E 

V^/CfAT 

Kinetic  energy/change  in 
thermal  energy 

Ekman 

Ek 

(r//2nL=')*/^  = 
(Ro/Re)*/^ 

(Viscous  force/Coriolis  force)^^^ 

Euler 

Eu 

Ap/pV'^ 

Pressure  drop  due  to  friction/ 
dynamic  pressure 

Froude 

Fr 

V/(sL)‘/2 

V/NL 

t(Inertial  force/gravitational  or 
buoyancy  force) 

Gay-Lussac 

Ga 

1/pAT 

Inverse  of  relative  change  in 
volume  during  heating 

Grashof 

Gr 

gh^fiATliP^ 

Buoyancy  force/viscous  force 

[Hall 

coefficient] 

Ch 

><lrL 

Gyrofrequency/ 
collision  frequency 

*(t)  Also  defined  as  the  inverse  (square)  of  the  quantity  shown. 
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Namo(s) 

Symbol 

Definition 

Significance 

Hartmann 

H 

= 

(RmRcC)*^^ 

(Magnetic  force/ 
dissipative  force) 

Kn\idscn 

Kn 

X/L 

Hydrodynamic  time/ 
collision  time 

Lewis 

Le 

kIV 

*Thermal  conduction/molecular 
diffusion 

Lorentz 

Lo 

V(c 

Magnitude  of  relativistic  effects 

Lundquist 

Lu 

f^oLVAfv  = 

A1  Rm 

J  X  B  force/resistive  magnetic 
diffusion  force 

Mach 

M 

V/Cs 

Magnitude  of  compressibility 
effects 

Magnetic 

Mach 

Mm 

1 

< 

II 

(Inertial  force/magnetic  force)^^^ 

Magnetic 

Reynolds 

Rm 

IM,LV/v 

Flow  velocity/magnetic  diffusion 
velocity 

Newton 

Nt 

F/pL^V^ 

Imposed  force/inertial  force 

Nussclt 

N 

aLfk 

Total  heat  transfer/thermal 
conduction 

Pdclet 

Pc 

LVfn 

Heat  convection/heat  conduction 

Poissenillc 

Po 

D^Ap/pLV 

Pressure  force/viscous  force 

Prandtl 

Pr 

uJk 

Momentum  diffusion/ 
heat  diffusion 

Rayleigh 

Ra 

gH^0AT/pK 

Buoyancy  force/diffusion  force 

Reynolds 

Re 

LVjv 

Inertial  force/viscous  force 

Richardson 

Ri 

(NHIAVf 

Buoyancy  effects/ 
vertical  shear  effects 

Rossby 

Ro 

V'/2nLsin  A 

Inertial  force/Coriolis  force 

Schmidt 

Sc 

ufV 

Momentum  diffusion/ 
molecular  diffusion 

Stanton 

St 

ocfpcpV 

Thermal  conduction  loss/ 
heat  capacity 

Stefan 

Sf 

aLT^/k 

Radiated  heat/conducted  heat 

Stokes 

S 

Viscous  damping  rate/ 
vibration  frequency 

Strouhal 

Sr 

Shiv 

Vibration  speed/flow  velocity 

Taylor 

Ta 

•(«/>') 

Centrifugal  force/viscous  force 
(Centrifugal  force/ 
viscous  force)^*^^ 

Timing, 

Boltzmann 

Th,  Bo 

pc„VleGT^ 

Convective  heat  transport/ 
radiative  heat  transport 

Weber 

W 

pLV^/T: 

Inertial  force/surface  tension 
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Nomenclature: 

B 

C,,c 

Cp 

D  =  2R 
F 

f 

9 

H,L 

k  =  pCpK 

JV= 

R 

T 

tl 

T 

V 

Va  =  B/(popy^^ 

Q 

r 

AH,  AV,  Ap,  AT 


V 

K,V 

A 

A 

fi~  pu 

MO 

u 

P 

p' 

S 

<7 

n 


Magnetic  induction 
Speeds  of  sound,  light 

Specific  heat  at  constant  pressure  (units  s“^  K“^) 

Pipe  diameter 

Imposed  force 

Vibration  frequency 

Gravitational  acceleration 

Vertical,  horizontal  length  scales 

Thermal  conductivity  (units  kgm“^s“^) 

Brunt-Vaisala  frequency 

Radius  of  pipe  or  channel 

Radius  of  curvature  of  pipe  or  channel 

Larmor  radius 

Temperature 

Characteristic  flow  velocity 
Alfv^n  speed 

dT 

Newton’s-law  heat  coefficient,  k-r—  =  aAT 

dx 

Volumetric  expansion  coefficient,  dVfV  —  0dT 

Bulk  modulus  (units  kgm“^s“^) 

Imposed  differences  in  two  radii,  velocities, 
pressures,  or  temperatures 

Surface  emissivity 
Electrical  resistivity 

Thermal,  molecular  diffusivities  (units  m^s“^) 
Latitude  of  point  on  earth’s  surface 
Collisional  mean  free  path 
Viscosity 

Permeability  of  free  space 
Kinematic  viscosity  (units  m^  s~^) 

Mass  density  of  fluid  medium 

Mass  density  of  bubble,  droplet,  or  moving  object 
Surface  tension  (units  kgs“^) 

Stefan-Boltzmann  constant 
Solid-body  rotational  angular  velocity 
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SHOCKS 

At  a  shock  front  propagating  in  a  magnetized  fluid  at  an  angle  6  with 
respect  to  the  magnetic  induction  B,  the  jump  conditions  are 

(1)  pU  =  pU  =  g; 

(2)  pU^+p-t-  Blllp  =  +  p  + 

(3)  pUV  —  B\\B±/li  =  pUV  —  B||Bx/p; 

(4)  B||  =  B||i 

(5)  C/Bx  -  VB||  =  BBx  -  VB||; 

(6)  ^(C/^  +  +  ui  +  -  VBjiBxj/ppW 

=  l(f/=  +  V*)  +  ®  +  (i/Sx^- V'B||Bx)/pp&. 

Here  U  and  V  are  components  of  the  fluid  velocity  normal  and  tangential  to 
the  front  in  the  shock  frame;  p  =  \/v  is  the  mass  density;  p  is  the  pressure; 
Bx  =  B  sin  9,  B||  =  B  cos  9;  p  is  the  magnetic  permeability  (p  =  47r  in  cgs 
units);  and  the  specific  enthalpy  is  in  =  c  +  pu,  where  the  specific  internal 

energy  e  satisfies  de  =  Tds  —  pdv  in  terms  of  the  temperature  T  and  the 

specific  entropy  a.  Quantities  in  the  region  behind  (downstream  from)  the 
front  are  distinguished  by  a  bar.  If  B  =  0,  then*® 

(7) U  -U  =  ((p-p){u-«)l'^’l 

(8)  (p-p)(u-0)-‘  =g’i 

(9)  ui  -  ti;  =  4 (p  -  p)(v  +  v); 

(10)  e  -  e  =  §(p  +  p)(v  -  t>). 

In  what  follows  wc  assume  that  the  fluid  is  a  perfect  gas  with  adiabatic  index 
^  =  1  +  2/n,  where  n  is  the  number  of  degrees  of  freedom.  Then  p  =  pHT/m, 
where  R  is  the  universal  gas  constant  and  m  is  the  molar  weight;  the  sound 
speed  is  given  by  =  {dpfdp),  =  ipv]  and  tu  =  76  =  7P^^/(7  -  1)- 
general  oblique  shock  in  a  perfect  gas  the  quantity  X  =  satisfies 

(11)  (X-l3/a)(X-cos**  ef  =  Xsin=9{ll  +  (r  -  \)/2cc]X  -  cos**  9},  where 
r  =  p/p,  Q  =  i  [7  +  1  -  (7  -  I)'-),  and  B  =  G.'^/Va^  =  in~ip/B^. 

The  density  ratio  is  bounded  by 

(12)  1  <  r  <  (7  +  1)/C7“  1)- 

If  the  shock  is  normal  to  B  (i.e.,  M  d  =  7r/2),  then 

(13)  =  (r/a){c,**  +  Vp**(l  +  (l-7/2)(r-  1)]}; 

(14)  U/U  =  B/B  =  r; 
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(15)  V  =  V-, 

(16)  p  =  p  +  (1  -  r-^)pU^  +  (1  -  r^)B^!2ti. 

If  0  =  0,  there  are  two  possibilities:  switch-on  shocks,  which  require  ^3  <  1  and 
for  which 

(17)  C/=  =:rVA^-, 

(18)  U  =  Va^/U-, 

(19)  =2B||^(r-l)(a-;S); 

(20)  V  =  UBj_IBif-, 

(21) p  =  p  +  pU^(l-a+p){l-T-'), 

and  acoustic  (hydrodynamic)  shocks,  for  which 

(22)  =  (r/a)Cs^\ 

(23)  U  =  U/v] 

(24)  V  =  B±  =0] 

(25)  p  =  p  +  pU^{l  -  r”'). 

For  acoustic  shocks  the  specific  volume  and  pressure  are  related  by 

(26)  v/v  =  [(7  +  1)P  +  (7  -  l)p|  /  ({7  -  1)P  +  (7  +  l)pl- 
In  terms  of  the  upstream  Mach  number  M  =  U/Cg, 

(27)  pIp  =  v/v  =  U/0  =  (7  +  1)M^/((7  -  1)M=  +  2]; 

(28)  p/p  =  (27 Af^  -  7  +  l)/(7  +  1); 

(29)  f/T  =  ((7  -  l)M^  +  2)(27M^  -  7  +  l)/(7  +  1)“M^ 

(30)  =  [(7  -  l)M^  +  2]/[27M2  -7  +  1]. 

The  entropy  change  across  the  shock  is 

(31)  As  =  s  -  s  =  c„  ln[(p/p)(p/p)'']. 

where  c^,  =  R/il  —  l)m  is  the  specific  heat  at  constant  volume;  here  R  is  the 
gas  constant.  In  the  weak-shock  limit  (M  —¥  1), 


(32)  As  ■ 


27(7  -  1) 

'  3(7+1) 


(M^ 


16'yR 
3(7  +  l)m 


{M-lf 


The  radius  at  time  t  of  a  strong  spherical  blast  wave  resulting  from  the  explo¬ 
sive  release  of  energy  in  a  medium  with  uniform  density  p  is 


(33)  Rs  =  Co(B^=/p)'''^ 

where  Co  is  a  constant  depending  on  7.  For  7  =  7/5,  Cq  =  1.033. 
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FUNDAMENTAL  PLASMA  PARAMETERS 

All  quantities  are  in  Gaussian  cgs  units  except  temperature  (T,  Te,  Ti) 
expressed  in  eV  and  ion  mass  (mi)  expressed  in  units  of  the  proton  mass, 
/i  =  mi/mp‘,  Z  is  charge  state;  k  is  Boltzmann’s  constant;  K  is  wavenumber; 
7  is  the  adiabatic  index;  In  A  is  the  Coulomb  logarithm. 

Frequencies 

electron  gyrofrcquency  fee  =  u»ce/2ir  =  2.80  x  10®BHz 

Wee  —  eB/rricC  =  1.76  x  lO^B  rad/sec 
ion  gyrofrcquency  fei  ~  a;ci/27r  =  1.52  X  10^Z/x“^BHz 

Wei  =  ZeB/rriiC  =  9.58  x  10®Z/j“^B  rad/sec 
electron  plasma  frequency  /pc  =  Wpc /Ztt  =  8.98  X  lO^rie  Hz 

Wpc  =  {A-nUe^  Irac)^''^ 

—  5.64  X  lO^rie^'^^  rad/sec 
ion  plasma  frequency  fpi  = 

=  2.10  X  Hz 

Wpi  =  (4jrniZ^e^/mi)^'^^ 

=  1.32  X  ]0®Z;i“*''^ni*'^^rad/sec 
electron  trapping  rate  t'Te  =  (cK’B/mc)^^^ 

=  7.26  X  sec"^ 

ion  trapping  rate  uri  =  {ZtKBIm.iY^'^ 

=  1.69  X  sec"^ 

electron  colli.sion  rate  Ve  ==  2.91  x  10“®ne  In  sec“^ 

ion  collision  rate  Vi  =4.80  x  In  ATj  sec  ^ 

Lengths 

electron  deBroglie  length  X  =  h/(mcA:Te)^^^  =  2.76  x  10“®re“^^^  cm 

classical  distance  of  e^fkT  =  1.44  x  10”^T“^  cm 

minimum  approach 

electron  gyroradius  r*  =  VTef<*fce  =  2.38Te^^^B  ^  cm 

ion  gyroradius  Vi  =  VTi/(^Jci 

=  1.02  X  cm 

plasma  skin  depth  cfwpe  =5.31  X  10®nc~^'^^cm 

Debye  length  Ap  =  (A:T/47rne^)^'^^ 

-7.43  X  10=^r^/^n"^^^cm 
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Velocities 

electron  thermal  velocity 

ion  thermal  velocity 

ion  sound  velocity 

Alfven  velocity 

Dimensionless 

(electron/proton  mass  ratio)*^^ 

number  of  particles  in 
Debye  sphere 

AIfv4n  velocity/speed  of  light 

electron  plasma/gyrofrequency 
ratio 

ion  plasma/gyrofrequency  ratio 
thermal/magnetic  energy  ratio 
magnetic/ion  rest  energy  ratio 

Miscellaneous 

Bohm  diffusion  coefficient 

transverse  Spitzer  resistivity 


=  4.19  X  10^ cm/sec 
VTi  =  (kTi/mi)^!'^ 

=  9.79  X  cm/sec 

C.  =  {■^ZkTJmiy^ 

=  9.79  X  10=(7ZTe/fi)‘/^  cm/sec 
VA  =  B  /  (4nnimi)^ 

=  2.18  X  cm/sec 

(me/rrip)^''^  =  2.33  X  10“^  =  1/42.9 
(4ir/3)nAD^  =  1.72  X 

=  3.21  X  lO-^ne'/^B-' 

u>p(/wc(  =  0.137^‘''2ni'/“B-' 

0  =  SnnkT/B^  =  4.03  X  lO-^nTB"* 

/STTniTTiiC^  = 

Db  =  (cfcT/16eB) 

=  6.25  X  10®rB-'  cmVsec 


V±  -  1.15  X  lO-^ZlnAT-^/'^sec 
=  1.03  X  10“'‘ZlnA7’“®/^ncm 
The  anomalous  collision  rate  due  to  low-frequency  ion-sound  turbulence  is 

!/•  X  Up,W/kT  =  5.64  X  lO‘^n,'^^W/kTsec~\ 


where  W  is  the  total  energy  of  waves  with  uj/K  <  VTt- 
Magnetic  pressure  is  given  by 

Bmag  =  B^/Stt  =  3,98  X  lo'’(B/So)^dynes/cm^  =  3,93(B/Bo)^  atm, 
where  Bq  =  10  kG  =  1  T. 

Detonation  energy  of  1  kiloton  of  high  explosive  is 

WkT  =  10‘^  cal  =  4,2  X  10'°  erg. 
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PLASMA  DISPERSION  FUNCTION 


Definition^®  (first  form  valid  only  for  Im^  >  0); 
_ _ _ 


r 

J  — o 


=  2iexp  J  dt  exp  ■ 


Physically,  (^  =  x  iy  is  the  ratio  of  wave  phase  velocity  to  thermal  velocity. 
Differential  equation: 

g  ^  -2  (1  +  CZ)  ,  Z(0)  =  g  +  2cg  +  2Z  =  0. 

Real  argument  (y  =  0): 

Z(x)  =  exp  (-1^)  -2  dt  exp  (t^)^  . 

Imaginary  argument  (x  =  0): 

Z(iy)  =  exp  (y^)  |1  -  erf(y)) . 

Power  series  (small  argument): 

ZiO  =  tTr'^^exp  (-<^)  -2<  (l  -2<^/3  +  4<Vi5  -  8CV105  +  •■•)  . 
Asymptotic  series,  |C|  )§>  1  (Ref.  17): 

ZiO  =  ix'/'^ffexp  (-C^)  -  <■'  (l  +  1/2C“  +  3/4C^  +  15/8C®  +  •••), 


f  0  y  >  kr' 

(  1  |y|  <  1*1  ‘ 

1.2  y  <  -|xr' 


Symmetry  properties  (the  asterisk  denotes  complex  conjugation): 

2(C*)  =  -1^(-C)]*; 

Z(C*)  =  [2(01  *  +  exp[-«*)*l  (y  >  0). 

Two-pole  approximations^®  (good  for  <  in  upper  half  plane  except  when  y  < 
exp(-x^),  X  1): 

0.50-h0.81i  0.50-0.81t  ^ 

Z(C)  - ,  a  =  0.51  —  0.81i; 

a-C  a*  +  C 

,  ,  0.50-t-0.96t  0.50-0.961  ^  . 


— — — ^ +  - - 6  =  0.48  -  0.91i. 

(6-C)2  ^  (fe*  +  C)" 
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COLLISIONS  AND  TRANSPORT 

Temperatures  are  in  eV;  the  corresponding  value  of  Boltzmann’s  constant 
is  A:  =  1.60  X  10“^^erg/eV;  masses  fx,  fx'  are  in  units  of  the  proton  mass; 
Ca  =  Za€  is  the  charge  of  species  a.  All  other  units  are  cgs  except  where 
noted. 

Relaxation  Rates 

Rates  are  associated  with  four  relaxation  processes  arising  from  the  in¬ 
teraction  of  test  particles  (labeled  a)  streaming  with  velocity  Vq  through  a 
background  of  field  particles  (labeled  j0): 


slowing  down 
transverse  diffusion 
parallel  diffusion 
energy  loss 


c,\0 

-  =  V/v 

dt 


where  the  averages  are  performed  over  an  ensemble  of  test  particles  and  a 
Maxwellian  field  particle  distribution.  The  exact  formulas  may  be  written^® 


=  2  [(1  -  1/“^'’; 


where 


=  47reQ^e;3^Ao^np/ma^Vo,^;  =  mpVa^  f2kT0\ 


dtp 
dx  ' 


and  Aq.^  =  InAc^  is  the  Coulomb  logarithm  (see  below).  Limiting  forms  of 
and  i/||  are  given  in  the  following  table.  All  the  expressions  shown 
have  units  cm®  sec“^.  Test  particle  energy  e  and  field  particle  temperature  T 
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are  both  in  eV;  /x  =  mi/mp  where  rrip  is  the  proton  mass;  Z  is  ion  charge 
state;  in  electron-electron  and  ion-ion  encounters,  field  particle  quantities  are 
distinguished  by  a  prime.  The  two  expressions  given  below  for  each  rate  hold 
for  very  slow  «  1)  and  very  fast  »  1)  test  particles,  respectively. 


Electron-electron 

uf'^'  «  5.8  X  10" 

«  5.8  X  10" 
«2.9  X  10- 

Electron-ion 

vl''' 

,/^'VniZ^Aei  «  2.5  X 
■.j'VniZ^Aci  K  1.2  X  10- 
Ion-electron 

«  1.6  X  10- 
i/^^'/n.Z^Ai,  K  3.2  X  10“ 
^'^'/neZ^Aio  w  1-6  X  10- 


•  7.7  X  lO-^e"^^^ 

.  7.7  X  10-^-^/" 

.  3.9  X  10-®re-'‘^^ 


.  1.7xl0-V^V^/= 


nj/Z^Z'^Aij/ 


Vm  mV 


Tij/Z^Z'^A,,/ 


Tl;/ Z^  Z^^  Aj o  1/5  / 1„  5/2 

■  — >  9.0  X  lO'V  '  M  Ti 

In  the  same  limits,  the  energy  transfer  rate  follows  from  the  identity 


except  for  the  case  of  fast  electrons  or  fast  ions  scattered  by  ions,  where  the 
leading  terms  cancel.  Then  the  appropriate  forms  are 


■  4.2  X  10  np 


_  8.9  X  10\,r/T)''V-‘  exp(-1836,i£/T)]  sec  ' 
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and 


— >  1.8  X 

-  1.1(//J')‘''’£“‘ exp(-fi'£/T)]  sec"‘. 


In  general,  the  energy  transfer  rate  is  positive  for  £  >  £q*  and  nega- 

tive  for  €  <  £«*,  where  x*  =  (m0\ma)ec,* /Tp  is  the  solution  of  -tli'{x*)  = 
(ma\m^)t/>(x*).  The  ratio  ^a*/1'0  is  given  for  a  number  of  specific  oc,  /3  in  the 
following  table: 


a\p 

i\e  e\e,  i\i 

e\p 

e\D 

e\T,  e\He^ 

eVHe"" 

Tp 

1.5  0.98 

4.8  X  10“^ 

2.6  X  10“® 

1.8  X  10“^ 

1.4  X  10'® 

When  both  species  are  near  Maxwellian,  with  Ti  ^  Tg,  there  are  just 
two  characteristic  collision  rates.  For  Z  =  1, 


=  2.9  X  10~®nATt"^'^sec“'i 
Ui  =4.8  X  sec"‘. 


Temperature  Isotropization 

Isotropization  is  described  by 


dTj.  _  1  dT||  _ 


where,  H  A  =  -  1  >  0, 


rriQ  l/2(AT||)3/2 


-3  +  (yl  +  3) 


tan-*(yl*/2) 

.41/2 


If  >1  <  0,  tan  is  replaced  by  tanh  ^  /{—A)^^^ .  For 

Tx  «  T||  =  r, 


!/J  =  8.2  X  10“''nAT“^‘'^sec"‘| 

=  1.9  X  10“®nAZ^/x"*^^T“^'^sec"'. 
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Thermal  Equilibration 

If  the  components  of  a  plasma  have  different  temperatures,  but  no  rela¬ 
tive  drift,  equilibration  is  described  by 


dTc 

dt 


-T^), 

/9 


where  ,  , 

=  1.8  X 

«  {mo.Tfi  +  mj3Tay/^ 

For  electrons  and  ions  with  Ti  =  T,  this  implies 

=  3.2  X  10-®Z^A/l»T^^^cm^  sec"*. 


Coulomb  Logarithm 

For  test  particles  of  mass  Tnc^  and  charge  Bq  =  ZqS  scattering  off  field 
particles  of  mass  mp  and  charge  ep  =  Zpe,  the  Coulomb  logarithm  is  defined 
as  A  =  In  A  =  ln(rnmx/rmiit)-  Here  rmin  is  the  larger  of  Baep/mapu^  and 
h/2mapu,  averaged  over  both  particle  velocity  distributions,  where  map  = 
mRm;3/(mo +mj3)  and  u  =  Vo —v^;  rmax  =  (^tt  /kT^)  where 

the  summation  extends  over  all  species  7  for  which  ~  kT-f/m-f.  If 

this  inequality  cannot  be  satisfied,  or  if  either  uwco  ^  <  ^max  or  uuicp  < 
Tnmxf  the  theory  breaks  down.  Typically  A  ss  10-20.  Corrections  to  the  trans¬ 
port  coefficients  are  0(A"M;  hence  the  theory  is  good  only  to  ~  10%  and  fails 

when  A  1.  .  .  ,  .  x 

The  following  cases  are  of  particular  interest: 

(a)  Thermal  electron-electron  collisions 


Ac.:  =  23  - 

=  24-ln(n/^=T.“'), 


r,  <  10  eV; 
Te  k,  10  eV. 


(b)  Electron-ion  collisions 

,  =23-ln(n//=Zn-=/"), 
=  24  -  In 
=  30 -In 


Aei  =  Ai, 
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Timjrrii  <  T,  <  lOZ’eV; 
Tivne/mi  <  lOZ^  eV  <  Te 
Te  <  TiZmdmi. 


(c)  Mixed  ion-ion  collisions 


X..,  =  X.,.  =  23  -  In 


mTj/  +  li'Ti  V  Ti  Til  ) 


(d)  Counterstreaming  ions  (relative  velocity  vd  =  ^dc)  in  the  presence  of 
warm  electrons,  kTi/rrii,  kT^i /m^/  <  vd^  <  kTefrUe 


Xii,  =  Xin 


35  -  In 


Fokker-Planck  Equation 

Dt  dt 


+  V  -  V/“  +  F  •  V. 


where  F  is  an  external  force  field.  The  general  form  of  the  collision  integral  is 

(a/“/s«)=oM  =  -  Efl 


■  J  —  uu)u”® 

■  I — r(v)v,,/^(v')  -  J-/^(v^)Vvr  (v)| 

L  tn^  me,  ) 


(Landau  form)  where  u  =  v'  —  v  and  /  is  the  unit  dyad,  or  alternatively, 

=  47rA„5  |r(v)VvH(v)  -  iVv  .  [/“(v)VvVvGCv)]  |  . 

where  the  Rosenbluth  potentials  are 


G(v)  =  J"  f^{v')ud%' 

H(v)  =  y /®(v')u-i<iV. 
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If  species  a  is  a  weak  beam  (number  and  energy  density  small  compared  with 
background)  streaming  through  a  Maxwellian  plasma,  then 


rriQ  +  m0 
1 
4 


v/“  -  : 


,  o\5 


2  « 

(k^/- vv)  •  Vv/“. 


VV  •  Vv/'^ 


B-G-K  Collision  Operator 

For  distribution  functions  with  no  large  gradients  in  velocity  space,  the 
Fokker-Planck  collision  terms  can  be  approximated  according  to 


Df. 

Dt 


=  u^AFc-fc)  +  i'ci(Fc-fc)-, 


Dfi 

Dt 


:Uic(Fi-fi)  +  ^ii(Fi-fi). 


The  respective  slowing-down  rates  i/f 


i  in  the  Relaxation  Rate  section 


above  can  be  used  for  1^0$,  assuming  slow  ions  and  fast  electrons,  with  e  re¬ 
placed  by  r«.  (For  Ucc  and  uu,  one  can  equally  well  use  and  the  result 
is  insensitive  to  whether  the  slow-  or  fast-test-particle  limit  is  employed.)  The 
Maxwellians  Fa  and  Fa  are  given  by 


Fc  =  "c.  ( 


Fc 


/  mo 

V  3/2 

exp  |- 

rme.(v  -  V„)^1  1 

V  2'^kTa , 

) 

L  2fcTe.  J  J 

/  mo 

,3/2 

exp 

rmc{v-v„)“11 

V  27rA:'ro  , 

) 

L  2fcf„  J  J 

where  na,  Va  and  Ta  are  the  number  density,  mean  drift  velocity,  and  effective 
temperature  obtained  by  taking  moments  of  fa-  Sorne  latitude  in  the  definition 
of  Ta  and  is  possible;^®  one  choice  is  Te  =  T*,  Ti  =  Te,  Ve  =  Vi,  Vj  =  Vg. 

Transport  Coefficients 

Transport  equations  for  a  multispecies  plasma: 


d^ng 

dt 


+  rio  V  •  Vo,  =  0; 


maTT-a  — 7“—  =  “ Vpa  —  V  ■  Pc 
dt 


+  Zgeug  E  +  -Vo  X  B  +  Rq; 

L  ^ 


36 


v„=-V-  q„-Po  ;  Vv„+Q„. 

2  at 

Here  d'^ /dt  =  d/dt  +  Va  •  V;  pa  =  nakTc,  where  k  is  Boltzmann’s  constant; 
Ra  =  “  Xly?  Qo0,  where  and  Qa0  are  respectively 

the  momentum  and  energy  gained  by  the  ath  species  through  collisions  with 
the  /3th;  Pq  is  the  stress  tensor;  and  qa  is  the  heat  flow. 

The  transport  coefficients  in  a  simple  two-component  plasma  (electrons 
and  singly  charged  ions)  are  tabulated  below.  Here  \\  and  ±  refer  to  the  di¬ 
rection  of  the  magnetic  field  B  =  bB;  u  =  —  v;  is  the  relative  streaming 
velocity;  rig  =  rii  =  n;  j  =  — neu  is  the  current;  uice  =  1.76  x  lO^Bsec”^  and 
ujci  =  {mc/mi)u}ce  are  the  electron  and  ion  gyrofrequencies,  respectively;  and 
the  basic  collisional  times  are  taken  to  be 


3.44X10°^ 

where  A  is  the  Coulomb  logarithm,  and 


sec, 


4>/7fn  Ae'* 


=  2.09  X  10^ 


J..3/2 

nA 


sec. 


In  the  limit  of  large  fields  (o/caTa  1,  a  =  i,  e)  the  transport  processes  may 
be  summarized  as  follows:^^ 


momentum  transfer 
frictional  force 

electrical 

conductivities 

thermal  force 


Rei=  — R-ic  =:  R  ==  Ru  d"  Rt> 

Ru  =  ne(j||/<7||  +jx/<^x); 

(T||  =  1.96crx;  =  ne^Tt/me] 

3u 

Rt=  -0.71nV||(fcre) - b  X  Vx(fcre); 

2uJceTe 


ion  heating 
electron  heating 


^  3mc  nk  ^  ^ 

Q.  = - i_(7;  _T.). 

rrxi  Tg 

-  R  .  u; 


ion  heat  flux 

ion  thermal 
conductivities 
electron  heat  flux 

frictional  heat  flux 


q*  =  — x(kTi)  d-  x  V ±{kTi)-, 

i  nkTiTi  i  2nkTi  ^  5nkTi 

«t|  =3.9 - ;  Kx  =  - 2 — !  =  n - ! 

*'  m*  2miWci 

Qe  =  Qu  d"  Qt*; 

,  ^  3nfcT'e 

q„  =  0.71nfcreU||  d - b  X  ux; 

SWccTe 
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thermal  gradient  =  — KjjV||(fc7V)  —  ±{kTe)  —  b  x  Vx(fcT'e); 

heat  flux 

e  nkTeTc  e  nkT^  e  5nkTe 

electron  thermal  «ll  =  3.2 - ;  Kj_  =  4.7 - - — ;  =  - 

conductivities  TncW^^rc  2meu;ce 

stress  tensor  (either  Pxx—  — —  +  ^w)  —  —  ’^yy)  — 

species)  2  2 

Pyy=  - yfW'-  +  ^yy)  +  -  ^yy'>  + 

p.„=  p,,  =  -mw,y  +  -  w'vv); 

P..=  P„  =  -7)2W'*.  -  7)4 Wv,; 

J°yz=  .PiW  =  4-  r]4Wxx-, 

(here  the  z  axis  is  defined  parallel  to  B); 

i  ^  3nkTi  i  6nkTi 

ion  viscosity  rj^  =  0.96nfcTjrj;  n,  =  - s — ;  ri,  =  - = — ; 

^  lOujjTi  5a»^fri 

*  -  i  _ 

^  2u>c»  "  lOci 


electron  viscosity  =0.7ZnkTeTc;  »?j=0.51- 

4  ^  nfcT*  e  _  nkTe 

2a»<;c  Wee 

For  both  species  the  rate-of-strain  tensor  is  defined  as 
dvi  dvk  2 

When  B  =  0  the  following  simplifications  occur: 


=2.0- 


Ru  =  nej/aii;  Rt  =  — 0.71nV(fcTe);  qi  =  — «m  V(/cTi); 


q*  =  0.71nfcTcU;  q^  =  — KjjV(fcT*e);  Pjk  =  —Vo^jk- 

For  u»ceTe  ^  1  WetTj,  the  electrons  obey  the  high-field  expressions  and  the 
ions  obey  the  zero-field  expressions. 

Collisional  transport  theory  is  applicable  when  (1)  macroscopic  time  rates 
of  change  satisfy  d/dt  ^  I/t,  where  r  is  the  longest  collisional  time  scale,  and 
(in  the  absence  of  a  magnetic  field)  (2)  macroscopic  length  scales  L  satisfy  L 
I,  where  1  =  ur  is  the  mean  free  path.  In  a  strong  field,  ojee'''  ^  I?  condition 
(2)  is  replaced  by  L^\  3>  I  and  L±  v/Ir*  {L±  re  in  a  uniform  field), 
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where  L\\  is  a  macroscopic  scale  parallel  to  the  field  B  and  L±_  is  the  smaller 
of  B/lVj.B|  and  the  transverse  plasma  dimension.  In  addition,  the  standard 
transport  coefficients  are  valid  only  when  (3)  the  Coulomb  logarithm  satisfies 
A  ^  1;  (4)  the  electron  gyroradius  satisfies  Ad,  or  STrneTTieC^  B^;  (5) 

relative  drifts  u  =  Vq  —  between  two  species  are  small  compared  with  the 
thermal  velocities,  i.e.,  kTa/rric,,  kT^Irct^-,  and  (6)  anomalous  transport 

processes  owing  to  microinstabilities  are  negligible. 

Weakly  Ionized  Plasmas 

Collision  frequency  for  scattering  of  charged  particles  of  species  a  by 
neutrals  is 

=  rioaf^°{kTa/maV^^ , 

where  no  is  the  neutral  density  and  is  the  cross  section,  typically 

5  X  10“^®  cm^  and  weakly  dependent  on  temperature. 

When  the  system  is  small  compared  with  a  Debye  length,  L  Ad,  the 
charged  particle  diffusion  coefficients  are 

Dci  =  kT^^  /irtcUoi , 

In  the  opposite  limit,  both  species  diffuse  at  the  ambipolar  rate 
^  —  f^e^i  _  (^t  "I" 

fix  “  fie  TiBg  +  TgDi 

where  fia  =  Ca/mat'o  is  the  mobility.  The  conductivity  a  a  satisfies  Co  = 

f^oi  • 

In  the  presence  of  a  magnetic  field  B  the  scalars  fi  and  ct  become  tensors, 
J“  =  O'"  ■  E  =  oj^B||  +  o“Ej,  +  o“E  X  b, 

where  b  =  B/B  and 

a  2  I 

^11  —  TIqCq 

ot  ct  ^  //  2,  2  \ 

<Tj_  =  <7,1  1^0  /(l/a 

O  Or  //  2  ,  2  \ 

=  ^11  Va'^cctlKya  +<^ce.)' 

Here  <tx  and  (Ta  are  the  Pedersen  and  Hall  conductivities,  respectively. 
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APPROXIMATE  MAGNITUDES 
IN  SOME  TYPICAL  PLASMAS 


Plasma.  Type 

n  cm  ^ 

T  eV 

Ufpe  sec  ^ 

\p  cm 

riXp^ 

Uei  sec  ^ 

Interstellar  gas 

1 

1 

6  X  10'* 

7  X  lO” 

4  X  10® 

7  X  10“” 

Gascons  nebula 

10^ 

1 

2  X  10° 

20 

lO” 

6  X  10“” 

Solar  Corona 

10” 

10^ 

2  X  10° 

2  X  10“' 

8  X  10° 

60 

DifTuse  hot  plasma 

10'” 

10^ 

6  X  10'° 

7  X  10“” 

4  X  10” 

40 

Solar  atmosphere, 

lo’"' 

1 

6  X  10" 

7  X  10“° 

40 

2  X  10° 

gas  discharge 

Warm  plasma 

10"* 

10 

6  X  10" 

2  X  10““ 

10” 

10^ 

Hot  plasma 

lO- 

10^ 

6  X  10" 

7  X  10'“ 

4  X  10“ 

4  X  10° 

Thermonuclear 

10l5 

lO'* 

2  X  10'” 

2  X  10'” 

10^ 

5  X  10“ 

plasma 

Theta  pinch 

10'“ 

10^ 

6  X  10'” 

7  X  10'” 

4  X  10” 

3  X  10® 

Dense  hot  plasma 

10'® 

10^ 

6  X  10'” 

7  X  10"° 

4  X  10” 

2  X  10'° 

Laser  Plasma 

10”” 

10^ 

6  X  10'“ 

7  X  10'” 

40 

2  X  10'” 

The  diagrjim  (faring)  gives  comparable  information  in  graphical  form. 
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50%  ionization 
hydrogen  plasma  / 


(4Tr/3)nX.  <  1 


1  2 
logio  T  (ev) 


-1 


0 


3 
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IONOSPHERIC  PARAMETERS'^ 


The  following  tables  give  average  nighttime  values.  Where  two  numbers 
arc  entered,  the  first  refers  to  the  lower  and  the  second  to  the  upper  portion 
of  the  layer. 


Quantity 

E  Region 

F  Region 

Altitude  (km) 

90-160 

160-500 

Number  density  (m~^) 

1.5  X  10'“-3.0  X  10'“ 

5  X  10'“-2  X  10" 

Height-integrated  number 
density  (m"^) 

9  X  lO*'’ 

4.5  X  10'“ 

Ion-neutral  collision 
frequency  (sec~^) 

2  X  10“-10“ 

0.5-0.05 

Ion  gyro-/collision 
frequency  ratio  Ki 

0.09-2.0 

4.6  X  10“-5.0  X  10“ 

Ion  Pederson  factor 

Ki/(1  + 

0.09-0.5 

2.2  X  10"“-2  X  I0-‘‘ 

Ion  Hall  factor 

KiV(l  +  «(“) 

8  X  10-'-0.8 

1.0 

Electron-neutral  collision 
frequency 

1.5  X  10‘'-9.0  X  10“ 

80-10 

Electron  gyro-/collision 
frequency  ratio 

4.1  X  10“-6.9  X  10“ 

7.8  X  10“ -6. 2  X  10“ 

Electron  Pcdcr.sen  factor 

Kc/(1  + 

2.7  X  10-“-1.5  X  10-’’ 

10-“-1.5  X  10-“ 

Electron  Hall  factor 

1.0 

1.0 

Mean  molecular  weight 

28-26 

22-16 

Ion  gyrofrcqiiency  (sec“^) 

180-190 

230-300 

Neutral  diffusion 

coefficient  (m^scc"^) 

30-5  X  10“ 

10“ 

The  terrestrial  magnetic  field  in  the  lower  ionosphere  at  equatorial  latti- 
tudes  is  approximately  Bq  —  0.35  x  10  ^  tesla.  The  earth’s  radius  is  Re  =  6371 
km, 
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SOLAR  PHYSICS  PARAMETERS^"* 


Parameter 

Symbol 

Value 

Units 

Total  mass 

Me 

1.99  X  10**** 

g 

Radius 

Rq  j 

6.96  X  10*“ 

cm 

Surface  gravity 

90 

2.74  X  10“ 

cm  s“^ 

Escape  speed 

■Doo 

6.18  X  10*' 

cms“^ 

Upward  mass  flux  in  spicules 

— 

1.6  X  10“® 

gcm“^  s“^ 

Vertically  integrated  atmospheric  density 

~ 

4.28 

gcm"^ 

Sunspot  magnetic  field  strength 

•^max 

2500-3500 

G 

Surface  effective  temperature 

To 

5770 

K 

Radiant  power 

Cq 

3.83  X  lO**** 

ergs“‘ 

Radiant  flux  density 

T 

6.28  X  10‘“ 

ergcm”^s“^ 

Optical  depth  at  500  nm,  measured 

T6 

0.99 

from  photosphere 

Astronomical  unit  (radius  of  earth’s  orbit) 

AU 

1.50  X  10*^ 

cm 

Solar  constant  (intensity  at  lAU) 

/ 

1.36  X  10® 

erg  cm“* 

Chromosphere  and  Corona^® 


Parameter  (Units) 

Quiet 

Sun 

Coronal 

Hole 

Active 

Region 

Chromospheric  radiation  losses 
(ergcm"^  s"^) 

Low  chromosphere 

2  X  10® 

2  X  10® 

^  10^ 

Middle  chromosphere 

2  X  10® 

2  X  10® 

lO"^ 

Upper  chromosphere 

3  X  10® 

3  X  10® 

2  X  10® 

Total 

4  X  10® 

4  X  10® 

^  2  X  10^ 

Transition  layer  pressure  (dyne  cm”®) 

0.2 

0.07 

2 

Coronal  temperature  (K)  at  1.1  Rq 

1. 1-1.6  X  10® 

10® 

2.5  X  10® 

Coronal  energy  losses  (erg cm”®  s”*) 

Conduction 

2  X  10® 

6  X  10“ 

10®-10^ 

Radiation 

10® 

10“ 

5  X  10® 

Solar  Wind 

^  5  X  10“ 

7  X  10® 

<  10® 

Total 

3  X  10® 

8  X  10® 

10^ 

Solar  wind  mass  loss  (gcm”®s”^) 

<  2  X  10“** 

2  X  10“*® 

<  4  X  10“** 
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THERMONUCLEAR  FUSION^" 


Natural  abundance  of  isotopes 
hydrogen 


Mass  ratios: 


npfn„  =  1.5  X  10"'* 
helium  ^  10"® 

lithium  TiLje/iiLjT  =  0.08 

mc/rriD  =  2.72  x  10“'*  =  1/3670 
1-65  x  10"^  =  1/60.6 
m^/rriT  =  1.82  x  lO""*  =  1/5496 
(mJrriTy^^^i  1.35  x  lO'^  =  1/74.1 


Absorbed  radiation  dose  is  measured  in  rads:  1  rad  =  10^  erg  g  The  curie 
(abbreviated  Ci)  is  a  measure  of  radioactivity:  1  curie  =  3.7x  10™  counts  sec"^ . 

Fusion  reactions  (branching  ratios  are  correct  for  energies  near  the  cross  section 
peaks;  a  negative  yield  means  the  reaction  is  endothermic):^^ 


(la) 

D  +  D 

-gg^T(1.01  MeV)  +  p(3.02  MeV) 

(lb) 

-=,:3B^He®(0.82McV)  +  n(2.45MeV) 
50% 

(2) 

D  +  T 

— — *He‘'(3.6McV)  +  n(14.1MeV) 

(3) 

D  -1-  He® 

- ►He''(3.6McV)  +  p{14.7MeV) 

(4) 

T  +  T 

- ►He''  +  2n  +  11,3  MeV 

(6a) 

He®  +  T 

-gy^Ho"  +  p  +  n  +  12.1  MeV 

(6b) 

.„„lHc'‘(4.8MeV)  +  D(9.6MeV) 
43% 

(6c) 

gy  lHc®(2.4MeV)  +  p(n.9MeV) 

(6) 

p  +  Li® 

- ►Hc‘‘(1.7MeV)  +  He®(2.3MeV) 

(7a) 

P  +  Li^ 

^2He^  +  17.3MeV 

(7b) 

-^Be'  +  n  1.6  MeV 

(8) 

D  +  Li® 

- ►2He'‘  +  22.4  MeV 

(9) 

P  +  B“ 

- >3  He''  +  8.7  MeV 

(10) 

n  +  Li® 

- *Hc''(2.1MeV)  +  T(2.7MeV) 

The  total  cross  section  in  barns  (1  barn  =  lO"^**  cm^)  as  a  function  of  E,  the 
energy  in  keV  of  the  incident  particle  [the  first  ion  on  the  left  side  of  Eqs. 
(l)-(5)],  assuming  the  target  ion  at  rest,  can  be  fitted  by^® 


ari^)  = 


As  +  [(A4  -  AiE)^  +  l]  ‘  A2 

E  [exp(AiB-i/2)  -  l] 


44 


where  the  Duane  coefficients  Aj  for  the  principle  fusion  reactions  are  as  follows: 


D-D 

(la) 

D-D 

(lb) 

D-T 

(2) 

D-He** 

(3) 

T-T 

(4) 

T-He® 

(5a-c) 

Ai 

46.097 

47.88 

45.95 

89.27 

38.39 

123.1 

A2 

372 

482 

50200 

25900 

448 

11250 

As 

4.36  X  IQ-* 

3.08  X  10-* 

1.368  X  10“^ 

3.98  X  10"® 

1.02  X  10"® 

0 

A4 

1.220 

1.177 

1.076 

1.297 

2.09 

0 

As 

0 

0 

409 

_ 

647 

0 

0 

Reaction  rates  av  (in  cm®  sec  averaged  over  Maxwellian  distributions: 


Temperature 

(keV) 

D-D 

(la-!-  lb) 

D-T 

(2) 

D-He® 

(3) 

T-T 

(4) 

T-He^ 

(5a-c) 

1.0 

1.5  X  10"®® 

5.5  X  10"®’ 

10"®® 

3.3  X  10"®® 

10-28 

2.0 

5.4  X  10"®* 

2.6  X  10"*® 

1.4  X  10"®® 

7.1  X  10"®* 

10"®® 

5.0 

1.8  X  10"*® 

1.3  X  10"*® 

6.7  X  10"®* 

1.4  X  10"*® 

2.1  X  10"®® 

10.0 

1.2  X  10"*® 

1.1  X  10"'® 

2.3  X  10"*® 

7.2  X  10"*® 

1.2  X  10"®° 

20.0 

5.2  X  10"*® 

4.2  X  10"*® 

3.8  X  10"'® 

2.5  X  10"*® 

2.6  X  10"*® 

50.0 

2.1  X  10"*® 

8.7  X  10"*® 

5.4  X  10"*® 

8.7  X  10"*® 

5.3  X  10"*® 

100.0 

4.5  X  10"'® 

8.5  X  10"*® 

1.6  X  10"'® 

1.9  X  10"'® 

2.7  X  10"*® 

200.0 

8.8  X  10-*® 

6.3  X  10"*® 

2.4  X  10"'® 

4.2  X  10"*® 

9.2  X  10"*® 

500.0 

1.8  X  10-*® 

3.7  X  10"*® 

2.3  X  10"*® 

8.4  X  10"*® 

2.9  X  10"*® 

1000.0 

2.2  X  10"*® 

2.7  X  10"'® 

1.8  X  10"'® 

8.0  X  10"*® 

5.2  X  10"*® 

For  low  energies  (T  ^  25  keV)  the  data  may  be  represented  by 

(5T)dd  =  2.33  X  sec"‘; 


(5U)dj.  =3.68  X  10"’^T"^^^exp(-19.94T“'^^)cm^  sec"*, 
where  T  is  measured  in  keV. 

The  power  density  released  in  the  form  of  charged  particles  is 

Pdd  —  3.3  X  10~^®7r£)^(  au)D£j  wattcm  ®  (including  the  subsequent 

D-T  reaction); 

Pdt  —  5.6  X  10“^®n£>nT(^)cT  wattcm”®; 

■^DHeS  =  2.9  X  10“^^n£,njjp3(^)£,He3 
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RELATIVISTIC  ELECTRON  BEAMS 

Here  7  =  (1  —  is  the  relativistic  scaling  factor;  quantities  in 

analytic  formulas  are  expressed  in  SI  or  cgs  units,  as  indicated;  in  numerical 
formulas,  I  is  in  amperes  (A),  B  is  in  gauss  (G),  electron  linear  density  N  is 
in  cm“^,  and  temperature,  voltage  and  energy  are  in  MeV;  Pz  =  Vz/c\  k  is 
Boltzmann’s  constant. 

Relativistic  electron  gyroradius: 

re  =  ^^(7^  -  1)*'’“  (cgs)  =  1.70  X  10^(7^  -  cm. 

eB 

Relativistic  electron  energy: 


W  =  mc^7  =  O.51I7  MeV. 


Bennett  pinch  condition: 

l‘‘  =  2N/c(Te  +  Ti)c“  (cgs)  =  3.20  x  10’‘'W(Te  +  Tt)  A“. 
Alfv6n-Lawson  limit: 

=  (mc^/e)^j7  (cgs)  =  {i-rrmc/ fxoe)Pz'r  (SI)  =  1.70  x  10'*/3a7A. 
The  ratio  of  net  current  to  Ia  is 


_L 

Ia 


V 

7’ 


Here  u  =  Nr^  is  the  Budker  number,  where  r*  =  e^/mc*  =  2.82  x  10  ^^cm 
is  the  classical  electron  radius.  Beam  electron  number  density  is 


where  J  is  the 
cm), 

and 


nt  =  2.08  X  10®J/3  ‘  cm 

current  density  in  A  cm”®.  For  a  uniform  beam  of  radius 

nt  =  6.63  X  10® cm”®, 

2re  v 
a  7” 


a  (in 
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Child’s  law:  (non-relativistic)  space-charge-limited  current  density  between 
parallel  plates  with  voltage  drop  V  (in  MV)  and  separation  d  (in  cm)  is 

J  =  2.34  X  A 


The  saturated  parapotential  current  (magnetically  self-limited  flow  along  equi- 
potentials  in  pinched  diodes  and  transmission  lines)  is^® 


/p  =  8.5  X  lO^Glln  [7+ (7^  -  A, 


where  G  is  a  geometrical  factor  depending  on  the  diode  structure: 

for  parallel  plane  cathode  and  anode 
of  width  w,  separation  d; 

for  cylinders  of  radii  Ri  (inner)  and  R2  (outer); 

for  conical  cathode  of  radius  Rc,  maximum 
separation  do  (at  r  =  /?c)  from  plane  anode. 

For  -+  0  (7  — >  1),  both  Ia  and  Ip  vanish. 

The  condition  for  a  longitudinal  magnetic  field  B*  to  suppress  filamentation 
in  a  beam  of  current  density  J  (in  Acm”^)  is 

B,  >  47^,(7J)‘'“  G. 

Voltage  registered  by  Rogowski  coil  of  minor  cross-sectional  area  A,  n  turns, 
major  radius  a,  inductance  L,  external  resistance  R  and  capacitance  C  (all  in 
SI): 


externally  integrated  V  —  (l/RC)(nA;/oI/2fra); 

self-integrating  V  =  (R/L)(nA^oI/2'^(^)  —  Rlln. 

X-ray  production,  target  with  average  atomic  number  Z  {V  ^  MeV): 

7]  =  x-ray  power/beam  power  =  7  X  10~^ZV. 

X-ray  dose  at  1  meter  generated  by  an  e-beam  depositing  total  charge  Q 
coulombs  while  V  >  0.84Vmax  >n  material  with  charge  state  Z\ 

D  =  rads. 
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BEAM  INSTABILITIES' 


30 


Namo 

Conditions 

Saturation  Mechanism 

Electron- 

electron 

Vd  >  i  =  1,2 

Electron  trapping  until 

Vej  -  Vd 

Buncman 

Ki  >  {M/m)'/^Vi, 

Vd  >  K 

Electron  trapping  until 

Ve  ~  Vd 

Beam-plasma 

Trapping  of  beam  electrons 

Weak  beam- 
plasma 

Vb  <  (np/nt)‘/®Vi, 

Quasilinear  or  nonlinear 
(mode  coupling) 

Beam-plasma 

(hot-electron) 

Ve  >  Vb  >  n 

Quasilinear  or  nonlinear 

Ion  acoustic 

r«  »  Ti,  Kd  »  c. 

Quasilinear,  ion  tail  form¬ 
ation,  nonlinear  scattering, 
or  resonance  broadening. 

Anisotropic 

temperature 

(hydro) 

T„j,  >  2T„|| 

Isotropi7/ation 

Ion  cyclotron 

Vd  >  2QVi  (for 

T,  a  Ti) 

Ion  heating 

Beam-cyclotron 

(hydro) 

Vd  >  C. 

Resonance  broadening 

Modified  two- 
stream  (hydro) 

Vd  <  (1+I3)'^^Va, 

Vd  >  c. 

Trapping 

Ion-ion  (equal 
beams) 

U  <2{1+P)''^Va 

Ion  trapping 

Ion-ion  (equal 
beams) 

U  <  2Cs 

Ion  trapping 

For  nomonclaturo,  sec  p.  50. 
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Parameters  of  Most  Unstable  Mode 


Name 

Wave 

Group 

Growth  Rate 

Frequency 

Number 

Velocity 

1 

Electron- 

0.9-f 

0 

electron 

2 

Vd 

Buneman 

0.4  I  —  1  lOe 

We 

Vd 

Beam-plasma 

We- 

We 

Vb 

Weak  beam- 
plasma 

2np  \vj 

We 

We 

n 

svj* 

Beam-plasma 

(hot-electron) 

(nj  Vi,“' 

Vi, 

n 

Ion  acoustic 

Wi 

^~D 

c. 

Anisotropic 

u>c  cos  9 

ry 

v»± 

temperature 

(hydro) 

Ion  cyclotron 

O.lQf 

1.20; 

ry 

Beam-cyclotron 

O.THe 

nOe 

o.ixy 

&  V.,; 

(hydro) 

i.r^ 

^  C3 

Modified  two- 

O.QOjy 

^Vd 

stream 

2 

Vd 

2 

(hydro) 

1.2^ 

Ion-ion  (equal 

0.40^ 

0 

0 

beams) 

Ion-ion  (equal 

0.4w,- 

0 

beams) 

For  nomenclature,  see  p.  50. 


In  the  preceding  tables,  subscripts  e,  i,  d,  b,  p  stand  for  “electron,”  “ion,” 
“drift,”  “beam,"  and  “plasma,”  respectively.  Thermal  velocities  are  denoted 
by  a  bar.  In  addition,  the  following  are  used: 


m 

electron  mass 

re,ri 

M 

ion  mass 

0 

V 

velocity 

T 

temperature 

Va 

Tie,  Tli 

number  density 

n 

harmonic  number 

Qh 

Cs  = 

ion  sound  speed 

U>e  ,  Wf 

plasma  frequency 

U 

Ad 

Debye  length 

gyroradius 

plasma/magnetic  energy 
density  ratio 
Alfv4n  speed 
gyrofrequency 
hybrid  gyrofrequency, 

relative  drift  velocity  of 
two  ion  species 


LASERS 


System  Parameters 

Efficiencies  and  power  levels  are  approximately  state-of-the-art  (1990).®^ 


Typo 

Wavelength 

Efficiency 

Power  levels  available  (W) 

(pm) 

Pulsed 

cw 

CO2 

10.6 

0.01-0.02 

(pulsed) 

>  2  X  10'® 

>  10® 

CO 

5 

0.4 

>  10® 

>  100 

Holmium 

2.06 

0.03t-0.1J 

>  10^ 

30 

Iodine 

1.315 

0.003 

>  10'® 

- 

Nd-glass, 

YAG 

1.06 

0.001-0.06t 

>o.u 

~  10^"  (ten- 
beam  system) 

1-10® 

♦Color  center 

1-4 

10-’ 

>  10® 

1 

♦Vibronic  (Ti 
Sapphire) 

0. 7-0.9 

>  0.1  X  lip 

10® 

1-5 

Ruby 

0.6943 

<  10-’ 

10*° 

1 

He-Ne 

0.6328 

10"* 

- 

1-50x10“® 

♦Argon  ion 

0.45-0.60 

10-’ 

5  X  lo" 

1-20 

♦OPO 

0. 4-9.0 

>  0.1  X  rip 

10® 

1-5 

N2 

0.3371 

0.001-0.05 

ID 

0 

*0 

- 

•Dye 

0.3-1. 1 

10-® 

>  10® 

140 

Kr-F 

0.26 

0.08 

>  10® 

500 

Xenon 

0.175 

0.02 

>  10® 

- 

♦Tunable  sources  flamp-driven  Jdiode-driven 

YAG  stands  for  Yttrium-Aluminum  Garnet  and  OPO  for  Optical  Parametric 
Oscillator;  rjp  is  pump  laser  efficiency. 
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Formulas 


An  e-m  wave  with  k  ||  B  has  an  index  of  refraction  given  by 
n±  =  [1  -  , 

where  ±  refers  to  the  helicity.  The  rate  of  change  of  polarization  angle  0  as  a 
function  of  displacement  s  (Faraday  rotation)  is  given  by 

de/ds  =  (A:/2)(n_  -  n+)  =  2.36  x  10‘‘NS/“^  cm-'. 

where  N  is  the  electron  number  density,  B  is  the  field  strength,  and  /  is  the 
wave  frequency,  all  in  cgs. 

The  quiver  velocity  of  an  electron  in  an  e-m  field  of  angular  frequency  w 
is 

vq  ~  =  25.67''^^Ao  cmsec”^ 

in  terms  of  the  laser  flux  I  =  with  I  in  watt/cm^,  laser  wavelength 

Aq  in  ^m.  The  ratio  of  quiver  energy  to  thermal  energy  is 

W^u/W^h  =  m^V0^/2kT  =  1.81  x  IQ-'^Ao^J/T, 

where  T  is  given  in  eV.  For  example,  if  /  =  lO'^  Wcm”^,  Ao  =  1  Mrri)  T  = 
2keV,  then  lVq„/lVth  «0.1. 

Pondermotive  force: 


where 

Nc  =  1.1  X  10^‘Ao-''cm“®. 

For  uniform  illumination  of  a  lens  with  /-number  F,  the  diameter  d  at 
focus  (85%  of  the  energy)  and  the  depth  of  focus  I  (distance  to  first  zero  in 
intensity)  are  given  by 

dK2A4F\e/eDL  and  I  R!  ±2F'‘ Xe/OoL- 

Here  6  is  the  beam  divergence  containing  85%  of  energy  and  6dl  is  the 
diffraction-limited  divergence: 

9dl  =  2.44A/6, 

where  b  is  the  aperture.  These  formulas  are  modified  for  nonuniform  (such  as 
Gaussian)  illumination  of  the  lens  or  for  pathological  laser  profiles. 


51 


ATOMIC  PHYSICS  AND  RADIATION 


Energies  and  temperatures  are  in  eV;  all  other  units  are  cgs  except  where 
noted.  Z  is  the  charge  state  (Z  =  0  refers  to  a  neutral  atom);  the  subscript  e 
labels  electrons.  N  refers  to  number  density,  n  to  principal  quantum  number. 
Asterisk  superscripts  on  level  population  densities  denote  local  thermodynamic 
equilibrium  (LTE)  values.  Thus  N„*  is  the  LTE  number  density  of  atoms  (or 
ions)  in  level  n. 

Characteristic  atomic  collision  cross  section; 

(1)  TTOo^  =  8.80  X  10“”  cm^. 


Binding  energy  of  outer  electron  in  level  labelled  by  quantum  numbers  n,  1: 


(2) 


££,(n.i)  =  - 


Z'^E^ 


whore  =  13.6eV  is  the  hydrogen  ionization  energy  and  Aj  =  0.751  ®, 
1  ^  5,  is  the  quantum  defect. 

Excitation  and  Decay 

Cross  section  (Bethe  approximation)  for  electron  excitation  by  dipole 
allowed  transition  m  n  (Refs.  32,  33): 


(3) 


=  2.36  X  10' 


.izfnm9(n,m)  ^ 

cAEnm 


where  fnm  is  the  oscillator  strength,  g{n,Tn)  is  the  Gaunt  factor,  e  is  the 
incident  electron  energy,  and  AEnm  —  En  —  Em- 

Electron  excitation  rate  averaged  over  Maxwellian  velocity  distribution,  Xmn 
=  Nciamnv)  (Refs.  34,  35): 


(4) 


AEn 


1/2 


where  (p(n,Tn))  denotes  the  thermal  averaged  Gaunt  factor  (generally  ~  1  for 
atoms,  ~  0.2  for  ions). 
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Rate  for  electron  collisiona!  deexcitation: 

(5)  Vnm  = 

Here  Nm* / Xn*  =  (5m/5n)exp(A£?nm/T’e)  is  the  Boltzmann  relation  for  level 
population  densities,  where  gn  is  the  statistical  weight  of  level  n. 

Rate  for  spontaneous  decay  n  — ^  m  (Einstein  A  coefficient)^^ 

(6)  A„m  =  4.3  X  10''(s„/s„)/„„(AB„„)^sec“’, 


Intensity  emitted  per  unit  volume  from  the  transition  n  — >•  m  in  an  optically 
thin  plasma; 

(7)  /„m  =  1.6  X  watt/cm^. 


Condition  for  steady  state  in  a  corona  model: 

(8)  WoJVe(tTo„u)  =  AfnAnO. 

where  the  ground  state  is  labelled  by  a  zero  subscript. 

Hence  for  a  transition  n  m  in  ions,  where  (5(n,0))  ss  0.2, 


(9) 


Inm  =5.1  X  10“ 


■,fnn,9oKNo  fAE„„ 


)%xp(- 


AB„oN 

T.  ) 


watt 

cm® 


Ionization  and  Recombination 

In  a  general  time-dependent  situation  the  number  density  of  the  charge 
state  Z  satisfies 

(10) 


Here  S{Z)  is  the  ionization  rate.  The  recombination  rate  cx(Z)  has  the  form 
a(Z)  =  ocr(Z)  -h  Nea3(Z),  where  ocr  and  03  are  the  radiative  and  three-body 
recombination  rates,  respectively. 


dN(Z} 


dt 


Ne  -  S{Z)N{Z)  -  a{Z)N{Z) 


+SiZ  -  1)N{Z  -  1)  +  a(Z  +  1)N{Z  +  1) 
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Classical  ionization  cross-section®®  for  any  atomic  shell  j 


(11)  CTi  =  6  X  10  cm^. 

Here  bj  is  the  number  of  shell  electrons;  Uj  is  the  binding  energy  of  the  ejected 
electron;  x  =  e/tlj,  where  e  is  the  incident  electron  energy;  and  is  a  universal 
function  with  a  minimum  value  ffmin  ~  0.2  at  x  ft;  4. 

Ionization  from  ion  ground  state,  averaged  over  Maxwellian  electron  distribu¬ 
tion,  for  0.02  <  T^/E^  g  100  (Ref.  35): 


(12)  S(Z)  =  10" 


(T,/E^)' 


(BZ)3/2(6.0-t-T,/B£) 


'(-#) 


cm®/sec, 


where  is  the  ionization  energy. 

Electron-ion  radiative  recombination  rate  (e+ N{Z)  N{Z 1)  +  hi') 
for  TclZ"^  ^  400 eV  (Ref.  37): 


(13)  a,(Z)=5.2xl0-“'Z^|f^^  j^O.43  +  i  ln(B2./^r) 

-l-0.469(£:2,/7’«)"’^^j  cmVsec. 

For  leV  <  Te/Z'^  <  15  eV,  this  becomes  approximately®^ 

(14)  ar{Z)  =  2.7  X 


Collisional  (three-body)  recombination  rate  for  singly  ionized  plasma:®® 

(15)  03  =  8.75  X  10"^^r.:"''  ®cm®/sec. 


Photoionization  cross  section  for  ions  in  level  n,  1  (short-wavelength  limit): 

(16)  o-ph(Ti,  0  =  1.64  X  cm®, 

where  K  is  the  wavenumber  in  Rydbergs  (1  Rydberg  =  1.0974  x  10®  cm  ^). 
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Ionization  Equilibrium  Models 

Saha  equilibrium:^® 


(17) 


NcN^*{Z) 

-  1) 


„Zrr  3/2 

n21  9\  Je 


cm 


-3 


where  is  the  statistical  weight  for  level  n  of  charge  state  Z  and  E^{n,l) 
is  the  ionization  energy  of  the  neutral  atom  initially  in  level  (n,/),  given  by 
Eq.  (2). 

In  a  steady  state  at  high  electron  density, 


N^N*{Z)  5(Z-1) 

N*{Z-\)  ~  Q3 


a  function  only  of  T. 

Conditions  for  LTE:^® 

(a)  CoIIisional  and  radiative  excitation  rates  for  a  level  n  must  satisfy 

(19)  YnmZlOAnm- 


(b)  Electron  density  must  satisfy 

(20)  Ne  ^  7  X 

Steady  state  condition  in  corona  model: 


(21) 


N{Z  -  1) 

N(Z)  ~S(Z-\)' 


Corona  model  is  applicable  if^® 

(22)  <N^<  cm“^, 

where  tj  is  the  ionization  time. 
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Radiation 

N.  B.  Etiergies  and  temperatures  are  in  eV;  all  other  quantities  are  in 
cgs  units  except  where  noted.  Z  is  the  charge  state  {Z  =  0  refers  to  a  neutral 
atom);  the  subscript  c  labels  electrons.  N  is  number  density. 

Average  radiative  decay  rate  of  a  state  with  principal  quantum  number  n  is 

(23)  A„  =  ^  Anm  =  1-6  X  10'° see. 

m<n 

Natural  linewidth  (AE  in  eV): 

(24)  AE  At  =  h  =  4.14  X  10“‘'^eVscc, 
where  At  is  the  lifetime  of  the  line. 

Doppler  width: 

(25)  AX/X  =  7.7  X  10’®(T/m)’'’, 

whore  is  the  mass  of  the  emitting  atom  or  ion  scaled  by  the  proton  mass. 
Optical  depth  for  a  Doppler-broadened  line:®^ 

(2C)  t  =  3,52  X  10"‘®/nmA(McV*.-T)‘'^NL  =  6,4  x  10-®X(fi/T)'^“tVL, 

where  /„m  is  the  absorption  oscillator  strength,  A  is  the  wavelength,  and  L  is 
the  physical  depth  of  the  plasma;  M,  N,  and  T  are  the  mass,  number  density, 
and  temperature  of  the  absorber;  fj.  is  M  divided  by  the  proton  mass.  Optically 
thin  moans  t  <  1. 

Resonance  absorption  cross  section  at  center  of  line: 

(27)  =  5.6  X  IQ-'^A^/AX  cm^. 

Wien  displacement  law  (wavelength  of  maximum  black-body  emission); 

(28)  =  2.50  X  cm. 

Radi.ntion  from  the  surface  of  a  black  body  at  temperature  T: 

(29)  W  =  1.03  X  10'*T''  watt/cm^. 
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Bremsstrahlung  from  hydrogen-like  plasma:^® 

(30)  Pbt  =  1.69  X  ^  [z^W(Z)]  watt/cm^ 

where  the  sum  is  over  all  ionization  states  Z. 

Bremsstrahlung  optical  depth:^^ 

(31)  T  =  5.0  X 

where  g  fa  1.2  is  an  average  Gaunt  factor  and  L  is  the  physical  path  length. 

Inverse  bremsstrahlung  absorption  coefficient'*^  for  radiation  of  angular  fre¬ 
quency  u): 

(32)  K  =  3.1  X  10-''Zne“lnAT-’/°w-^(l  -  cm“'; 

here  A  is  the  electron  thermal  velocity  divided  by  V,  where  V  is  the  larger  of 
uj  and  Up  multiplied  by  the  larger  of  Ze^ fkT  and  h/{mkT)^^'^ . 

Recombination  (free-bound)  radiation: 

(33)  Pr  =  1.69  X  ^  j^Z^N{Z)  (%— )]  wa‘t/cm^ 


Cyclotron  radiation^®  in  magnetic  field  B: 

(34)  Pc  =  6.21  X  10"^®B^AreT«  watt/cm^ 

For  NgkTe  =  NikTi  =  B^/lQn  {0  ss  l,  isothermal  plasma),^® 

(35)  Pc  =  5.00  X  watt/cm®. 


radiation 


energy  loss  e-folding  time  for  a  single  electron:**^ 
9.0  X  10® 

tc  »  - r - sec, 


where  7  is  the  kinetic  plus  rest  energy  divided  by  the  rest  energy  mc^. 
Number  of  cyclotron  harmonics^*^  trapped  in  a  medium  of  finite  depth  L: 
(37)  =  (57^BL)*'®, 

where  /3  =  BivNkT/B^. 

Line  radiation  is  given  by  summing  Eq.  (9)  over  all  species  in  the  plasma. 
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ATOMIC  SPECTROSCOPY 


Spectroscopic  notation  combines  observational  and  theoretical  elements. 
Obscrvationally,  spectral  lines  are  grouped  in  scries  with  line  spacings  which 
decrease  toward  the  series  limit.  Every  line  can  be  related  theoretically  to  a 
transition  between  two  atomic  states,  each  identified  by  its  quantum  numbers. 

Ionization  levels  arc  indicated  by  roman  numerals.  Thus  C  I  is  unionized 
carbon,  C  II  is  singly  ionized,  etc.  The  state  of  a  one-electron  atom  (hydrogen) 
or  ion  (Hell,  Li  III,  etc.)  is  specified  by  identifying  the  principal  quantum 
number  n  =  1,2,...,  the  orbital  angular  momentum  /  =  0,  1, . . . ,  n  -  1,  and 
the  spin  angular  momentum  s  =  ±^.  The  total  angular  momentum  j  is  the 
magnitude  of  the  vector  sum  of  I  and  s,  j  =  I  ±  5  (j  >  J)-  The  letters  s, 
p,  d,  f,  g,  h,  i,  k,  1,  ...,  respectively,  are  associated  with  angular  momenta 
I  ’=  6,  1,  2,  3,  4,  5,  6,  7,  8,  ... .  The  atomic  states  of  hydrogen  and  hydrogenic 
ions  arc  degenerate:  neglecting  fine  structure,  their  energies  depend  only  on  n 
according  to 

RooheZ^n-^  _  KyZ'^ 

1  +  m/M 

whore  h  is  Planck's  constant,  c  is  the  velocity  of  light,  m  is  the  electron  mass, 
M  and  Z  arc  the  mass  and  charge  state  of  the  nucleus,  and 

Hoc  =  109,737  cm”* 

is  the  Rydberg  constant.  If  E„  is  divided  by  he,  the  result  is  in  wavenumber 
units.  The  energy  as.sociatcd  with  a  transition  m  — ^  n  is  given  by 

AiS„„  =  Ry{l/m^  -  1/n"), 

with  m  <  n  {m>  n)  for  absorption  (emission)  lines. 

For  hydrogen  and  hydrogenic  ions  the  series  of  lines  belonging  to  the 
transitions  m  -¥  n  have  conventional  names: 


Transition 

1  n 

2  n 

3  ^  n 

4  n 

5  — ►  n 

6  — >  n 

Name 

Lyman 

B  aimer 

Paschen 

Brackett 

Pfund 

Humphreys 

Successive  lines  in  any  series  are  denoted  a,  7,  etc.  Thus  the  transition  1 
3  gives  rise  to  the  Lyman-/?  line.  Relativistic  effects,  quantum  electrodynamic 
effects  (e.g.,  the  Lamb  shift),  and  interactions  between  the  nuclear  magnetic 
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moment  and  the  magnetic  field  due  to  the  electron  produce  small  shifts  and 
splittings,  ^  10“^cm“^;  these  last  are  called  “hyperfine  structure.” 

In  many-electron  atoms  the  electrons  are  grouped  in  closed  and  open 
shells,  with  spectroscopic  properties  determined  mainly  by  the  outer  shell. 
Shell  energies  depend  primarily  on  n;  the  shells  corresponding  to  n  =  1,  2, 
3, . . .  are  called  K,  L,  M,  etc.  A  shell  is  made  up  of  subshells  of  different 
angular  momenta,  each  labeled  according  to  the  values  of  n,  I,  and  the  number 
of  electrons  it  contains  out  of  the  maximum  possible  number,  2(21  +  1).  For 
example,  2p®  indicates  that  there  are  5  electrons  in  the  subshell  corresponding 
to  /  =  1  (denoted  by  p)  and  n  =  2. 

In  the  lighter  elements  the  electrons  fill  up  subshells  within  each  shell 
in  the  order  s,  p,  d,  etc.,  and  no  shell  acquires  electrons  until  the  lower  shells 
are  full.  In  the  heavier  elements  this  rule  does  not  always  hold.  But  if  a 
particular  subshell  is  filled  in  a  noble  gas,  then  the  same  subshell  is  filled  in 
the  atoms  of  all  elements  that  come  later  in  the  periodic  table.  The  ground 
state  configurations  of  the  noble  gases  are  as  follows: 

He  ls“ 

Ne  ls^2s^2p® 

At  ls^2s^2p®3s^3p® 

Kr  ls^2s^2p®3s^3p®3d‘'’4s^4p® 

Xe  ls^2s^2p«3s“3p®3d'‘’4s^4p®4d’°5s^6p® 

Rn  ls^2s*2p®3s^3p*3d'°4s^4p®4d’‘’4f'^5s^5p®5d'“6s'“6p® 

Alkali  metals  (Li,  Na,  K,  etc.)  resemble  hydrogen;  their  transitions  are  de¬ 

scribed  by  giving  n  and  I  in  the  initial  and  final  states  for  the  single  outer 
(valence)  electron. 

For  general  transitions  in  most  atoms  the  atomic  states  are  specified  in 
terms  of  the  parity  (—1)^^*  and  the  magnitudes  of  the  orbital  angular  momen¬ 
tum  L  =  Slj,  the  spin  S  =  Esj,  and  the  total  angular  momentum  J  =  L  -f  S, 
where  all  sums  are  carried  out  over  the  unfilled  subshells  (the  filled  ones  sum 
to  zero).  If  a  magnetic  field  is  present  the  projections  M^,  Ms,  and  M  of 
L,  S,  and  J  along  the  field  are  also  needed.  The  quantum  numbers  satisfy 
\Ml\  <  |Afs|  <  S  <  Z//2,  and  |Mj  <  J  <  L  +  S,  where  i'  is  the 

number  of  electrons  in  the  unfilled  subshell.  Upper-case  letters  S,  P,  D,  etc., 
stand  for  L  =  0,  1,  2,  etc.,  in  analogy  with  the  notation  for  a  single  electron, 
For  example,  the  ground  state  of  Cl  is  described  by  3p®  ^^3/2’  part 

indicates  that  there  are  5  electrons  in  the  subshell  corresponding  to  n  =  3  and 
/  =  1.  (The  closed  inner  subshells  ls^2s*2p®3s^,  identical  with  the  configura¬ 
tion  of  Mg,  are  usually  omitted.)  The  symbol  ‘P’  indicates  that  the  angular 
momenta  of  the  outer  electrons  combine  to  give  L  =  1.  The  prefix  ‘2’  repre¬ 
sents  the  value  of  the  multiplicity  25-1-1  (the  number  of  states  with  nearly  the 
same  energy),  which  is  equivalent  to  specifying  5  =  5.  The  subscript  3/2  is 
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the  value  of  J.  The  superscript  ‘o’  indicates  that  the  state  has  odd  parity;  it 
would  be  omitted  if  the  state  were  even. 

The  notation  for  excited  states  is  similar.  For  example,  helium  has  a  state 
ls2s  ^Si  which  lies  19.72eV  (159,856cm~^)  above  the  ground  state  Is^  ^Sq- 
But  the  two  “terms”  do  not  “combine”  (transitions  between  them  do  not  occur) 
because  this  would  violate,  e.g.,  the  quantum-mechanical  selection  rule  that 
the  parity  must  change  from  odd  to  even  or  from  even  to  odd.  For  electric 
dipole  transitions  (the  only  ones  possible  in  the  long-wavelength  limit),  other 
selection  rules  are  that  the  value  of  I  of  only  one  electron  can  change,  and  only 
by  M  =  ±1;  AS"  =  0;  AL  =  ±1  or  0;  and  AJ  =  ±1  or  0  (but  L  =  0  does  not 
combine  with  L  =  0  and  7  =  0  does  not  combine  with  j  =  0).  Transitions 
are  possible  between  the  helium  ground  state  (which  has  S'  =  0,  L  =  0,  J  =  0, 
and  even  parity)  and,  e.g.,  the  state  ls2p  (with  S  =  0,  L  =  1,  J  =  1, 
odd  parity,  excitation  energy  21.22  eV).  These  rules  hold  accurately  only  for 
light  atoms  in  the  absence  of  strong  electric  or  magnetic  fields.  Transitions 
that  obey  the  selection  rules  are  called  “allowed”;  those  that  do  not  are  called 
“forbidden." 

The  amount  of  information  needed  to  adequately  characterize  a  state  in¬ 
creases  with  the  number  of  electrons;  this  is  reflected  in  the  notation.  Thus^"* 
on  has  an  allowed  transition  between  the  states  2p^3p^ 
*F®y2  2p^(^D)3d'  ^F7/2  between  the  states  obtained  by  changing 

J  from  7/2  to  5/2  in  either  or  both  terms).  Here  both  states  have  two  elec¬ 
trons  with  n  =  2  and  I  =  1;  the  closed  subshells  ls^2s®  are  not  shown.  The 
outer  (n  =  3)  electron  has  i  =  1  in  the  first  state  and  /  =  2  in  the  second. 
The  prime  Indicates  that  if  the  outermost  electron  were  removed  by  ionization, 
the  resulting  ion  would  not  be  in  its  lowest  energy  state.  The  expression  (^D) 
give  the  multiplicity  and  total  angular  momentum  of  the  “parent”  term,  i.e., 
the  subshell  immediately  below  the  valence  subshell;  this  is  understood  to  be 
the  same  in  both  states.  (Grandparents,  etc.,  sometimes  have  to  be  specified 
in  heavier  atoms  and  ions.)  Another  example^^  is  the  allowed  transition  from 
2p^(^P)3p  ^Pj/2  ^^3/2)  2p^(*D)3<i^  ^Si/2,  in  which  there  is  a  “spin 

flip”  (from  antiparallel  to  parallel)  in  the  n  =  2,/  =  1  subshell,  as  well  as 
changes  from  one  state  to  the  other  in  the  value  of  I  for  the  valence  electron 
and  in  L. 

The  description  of  fine  structure,  Stark  and  Zeeman  effects,  spectra  of 
highly  ionized  or  heavy  atoms,  etc.,  is  more  complicated.  The  most  important 
difference  between  optical  and  X-ray  spectra  is  that  the  latter  involve  energy 
changes  of  the  inner  electrons  rather  than  the  outer  ones;  often  several  electrons 
participate. 
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AFTERWORD 
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Formulary  has  been  set  in  T^jX  by  Dave  Book,  Todd  Brun,  and  Robert  Scott. 
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